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Raytheon  Corporation,  and  the  FAA  Technical  Center  in  Atlantic  City,  NJ. 

The  staff  members  of  the  Weather  Sensing  Group  have  documented  their  studies  in 
four  major  areas;  Operational  Systems  (TDWR  Operational  Test  and  Evaluation  re¬ 
sults);  Radar  Operations  (future  airport  weather  surveillance  requirements,  a  “machine 
intelligent”  gust  front  detection  algorithm,  microburst  asymmetry  study  results,  a  shear  - 
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PRELIMINARY  RESULTS  OF  THE  WEATHER  TESTING  CXDMPONENT  OF  THE  TERMINAL 
DOPPLER  WEATHER  RADAR  OPERATIONAL  TEST  AND  EVALUATON 


Steven  V.  Vasiloff',  Michael  D.  Eilts',  Elizabeth  TuFcich^ 
Jim  Wiele/’,  Maik  Isanungei^,  and  Paul  Biron^ 

’National  Severe  Stonns  Laboratory 
Nonnan,  Oklahoma  73069 

federal  Aviation  Administration  Technical  Center 
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Atlantic  City,  New  Jersey  08405 

’Raytheon  Coq). 

Sudbury,  I^sachusetts  01776 

^Massachusetts  Institute  of  Technology/LiiKXtln  Laboratoty 
Lexington,  Massachusetts  02173 


1.  INTRODUCTION 

The  Terminal  Doppler  Weather  Radar  (TDWR) 
system  (Turnbull  et  al.  1989),  which  has  been 
developed  by  Raytheon  Co.  for  the  Federal  Aviation 
Administration  (FAA),  provides  automatic  detection 
of  microbursts  (Fujita  1985)  and  low-altitude  wind 
shear.  Microburst- and  gust  fiont-induced  wind  shear 
can  result  in  a  sudden,  large  change  in  airqreed  which 
can  have  a  disastrous  effect  on  turoaft  p^ormance 
during  take  off  or  landing.  The  snxxtd  major 
function  TDWR  b  to  improve  air  traffic 
management  through  forecasts  of  wind  shifts, 
precipitation  and  other  weather  hazards.  The  TDWR 
system  generates  Doppler  velocity,  reflectivity,  and 
spectrum  width  data.  The  base  daa  are  autoouaically 
dealiased  and  clutter  is  removed  through  filtering  and 
mapping.  Precipitationandwindshear  products,  such 
as  microbursts  and  gust  fronts,  are  displayed  as 
graphic  products  on  the  Geographic  Situation  Display 
(Fig.  1)  which  is  intended  for  use  by  Air  Tr^c 
Control  supervisors.  Alphanumeric  messages 
indicating  the  various  windshear  alerts  and  derived 
airspeed  losses  and  gains  are  sent  to  a  flat  panel 
ribt^  diqilay  which  is  by  die  controUos  in  the 
control  tower. 

The  TDWR  proof-of-conceptand  operational 
feasibility  have  been  demonstrated  in  a  number  of 
FAA-sponsored  tests  and  evaluadons  conducted  by 


Figure  1.  Black  and  white  reproduc&m  of  a 
geographic  situation  (Bsplay.  The  figure  it  centered 
near  Wiley  Post  Airport.  Precipitation  cm  be  seen 
to  the  west  cf  the  airport  and  is  grey  scaled  (darkest 
is  heaviest).  The  heavy  solid  line  is  the  gust  front 
and  the  dashed  lines  are  10  and  20-nun  forecasts.  A 
control  panel  to  the  right  is  not  shown  due  to  space 
constraints.  Range  rings  are  in  nautical  ntHss. 
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Massachusetts  Instuute  of  Technology’s 
Lincoln  Laboratory  (MIT/LL)  in  Memphis.  TN 
(1985);  HuntsviUe.  AL  (1986);  Denver,  CO  (1987, 
1988);  Kansas  City,  MO  (1989),  and  Orlando  FI. 
(1990-1992). 

In  Older  to  verify  that  the  TDWR  meets  FAA 
operational  suitability  and  effectiveness  requirements, 
an  Operational  Test  &  Evaluation  (OT&E)  was 
conducted  at  the  Oklahoma  City  site  during  the  period 
from  24  August  to  30  October  1992.  The  testing 
addressed  National  Ak^ce  System  (NAS)-SS-1000 
requirements,  weather  detection  performance,  safety, 
operational  system  performance,  maintenaiKe 
instruction  books.  Remote  Maintenance  Monitoring 
System  (RMMS),  system  adaptable  parameters, 
bullgear  wear,  and  limited  Air  Traffic  (AT)  suitability. 

The  TD^Tl  OT&E  Integration  and 
Operational  testing  was  conducted  using  a  variety  of 
methods  dependent  on  the  area  being  tested.  This 
paper  discusses  primarily  the  weather  detection 
performance  testing.  A  rough  analysis  was  performed 
on  the  algorithm  output  and  the  base  data  to 
determine  the  performance  of  the  TDWR  in  detecting 
wind  shear  phenomena.  Final  results  will  be  available 
after  additional  testing,  which  is  scheduled  for  Spring 
of  1993,  and  post  analysis  is  xmducied. 

2.  SYSTFM  DESCRIPTION 

The  TDWR  is  a  C-Band  weaflier  radar 
consisting  of  a  25-fooi  diameta,  center-fed  parabolic 
reflector  antenna,  with  the  feed  mounted  on  a  tripod. 
The  antenna  bcamwidtfi  is  0.55  degrees.  The  TDWR 
transmitter  tube  is  a  250  kW  peak-power  pulsed 
klystron,  transmitting  a  1.1  microsecond  pulse  (-6  dB 
width)  at  pulse  iq)idtion  frequencies  from  250  to 
2000  IIz.  Large  dynamic  range  is  provided  for  both 
good  clutter  suppression  and  accurate  reflectivity 
measurements.  The  total  dynamic  range,  over  128 
dB,  is  obtained  with  26  dB  of  Senritivity  Time 
Control  (STC),  42  dB  of  Automatic  Gm  Control 
(AGC),  and  circuits  linear  over  61  dB,  For  more 
detail  on  system  characteristics  the  reader  is  referred 
to  Wieler  arid  Shrader  1991,  and  Michelson  et  al. 
1990. 

3.  WEATHER  TESTING  RESULTS 

This  section  describes  results  of  the 
weather-detection  component  of  the  TDWR  OT&E. 
The  National  Severe  Storms  Laboratory  developed 
procedures  to  address  three  main  components  of  the 
TDV/R  weather  detecting  capalality:  1)  Assess  tftat 
the  base  data  are  of  Mgh  quality;  2)  Determine  if  the 


algorithms  are  functioning  properly;  3)  Verify  that 
appropriate  alarms  are  dissenunated  according  to  the 
system  design. 

3.1  Base  data 

In  a  qualitative  assessment,  the  judgement  of 
the  Investigative  Panel,  a  group  of  expert  radar 
meteorologists  from  the  National  Cen\CT  for 
Atmo^thcric  Research  (NCAR),  MITyLincoln 
Laboratory  and  NSSL,  is  that  the  raw  base  data  were 
of  high  quality.  All  the  NAS  requirements  related 
to  measurements  of  base  velocity,  reflectivity,  and 
spectrum  width  were  fulfilled. 

Suppression  of  ground  clutter  is  important  for 
"clean”  base  data  and  for  prevention  of  fate  alarms. 
The  TDWR  uses  two  techniques  to  reduce 
contamination  from  ground  clutter  sources.  The  first 
involves  the  use  of  IIR  (infinite  impute  response) 
Alters  to  suppress  high  reflectivity  returns  that  have 
near-zero  vdocities.  The  second  st^  is  to  generate  a 
clutter  residue  editing  map  (CREKQ  on  a  dear  day 
witfi  no  weather  echoes.  Experience  has  shown  the 
best  condition  for  taking  clutter  measurements  is  low 
clear-air  reflectivities  and  a  minimal  amount  of 
moving  sources  such  as  birds  aitd  insects. 

Clutter  map  generation  begins  with  an 
automated  determination  of  a  clear-air  reflectivity 
(CAR)  value.  The  CAR  estimate  and  a  maximum 
vekicify  threshold  are  then  used  to  distinguish  clutter 
residue  from  clear-tur  retunts  caused  by  birds  andAtr 
insects.  It  was  determined  that  the  process  for 
brulding  the  CREMs  may  not  be  straightforward  at 
each  site.  A  large  amount  of  bird  activity  (and 
posribly  insects)  as  well  as  evolving  boundary  layer 
characteristics  resulted  in  having  to  make  die  CAR 
estimate  before  sunrise.  It  is  anticipated  that  each  site 
will  have  its  own  peculiarities. 

32  Algorithm  performance 

32.1  Microburst  detection  algorithm 

The  Microburst  Detection  Algorithm  (MDA) 
detects  low-altitude  divergent  shears  associated  with 
storms  (Merritt  et  aL  1989).  Since  most  of  the  «toia 
were  collected  in  high-wind  environments,  classic 
microburst  signamres  were  rarely  seen. 

Because  a  limited  amount  of  data  are 
available  for  evaluation,  only  rough  assessments  of 
the  MDA  performance  are  available.  Detection 
peifonnance  was  excellent  although  a  number  of 
alarms  were  observed.  Large  flocks  of  birds  dqraiting 
from  roosting  sites  in  die  early  morning  hours  caused 
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divergent  signatures  similar  to  miciobursts.  There 
were  numerous  £fdse  aJarms  from  this  phenomena 
prior  to  the  implementation  of  the  storm  celt  test 
which  validates  alarms  based  on  reflectivity  aloft 
(Evans  1990X  Many  of  the  other  false  alarms  were 
the  result  of  clutter  breakthrough  or  noise  in  the 
velocity  data.  Note  the  clutter  breakthrough  around 
the  airport  in  Fig.  1. 

32.2  Gust  front  detection  algorithm 

The  Gust  Front  Detection  Algorithm  (GFDA) 
detects  lines  of  convergence  in  Doppler  radial  velocity 
fields  and  forecasts  the  movement  of  these  wind  shifts 
and  the  winds  betiind  them  up  to  20  minutes  in 
advance  (Hermes  et  aL  1992).  An  example  of  a  gust 
front  detection  is  shown  in  Fig.  1.  Test  results 
indicate  that  the  forecasting  function  of  the  algorithm 
performed  well,  as  did  the  estimation  of  winds  behind 
the  wind  shift.  As  a  gust  front  passed  over  Wiley 
Post  Airport  on  10  September,  data  from  the  local 
automatic  surface  obser^g  site  agreed  well  with  the 
windshift  value  behind  the  gust  front  detection. 

The  GFDA  performed  well  during  the 
OT&E  However,  as  with  the  MDA,  false  alarms 
were  observed.  The  majority  of  the  false  alarms  were 
detected  outside  20  km  and  thus  would  not  impact  the 
airport  Many  of  the  GFDA  false  alanns  were  the 
result  of  veni^  wind  shear  where  the  winds  veered 
(turned)  vrith  height  and  were  not  the  typical 
low-level  jet  atuations  which  have  been  noted  as  a 
potential  cause  of  false  alanns.  However,  some  false 
alarms  were  caused  by  more  classical  low-level  jets. 
We  are  classifying  these  vertical  wind  shears  and 
low-level  jets  as  false  alarms  since  they  were  not 
generated  by  a  thunderstorm  outflow.  It  could  be 
argued,  however,  that  these  features  are  operationally 
significanL  There  were  also  a  number  of  false  alanns 
due  to  dealiaring  errors  and  range  folding.  Minor 
software  changes  are  being  made  to  reduce  the 
number  of  false  alanns. 

4.  CONCLUSIONS 

The  initial  Operational  Test  &  Evaluation  of 
the  Temunal  Doppler  Radar  showed  that  weather 
detection  perftxmance  was  generally  acceptable.  A 
set  ofNAS-SS-KXX)  requirements  was  satisfied.  Base 
data  quality  appears  to  be  excellent  and  the  two 
primary  algorithms,  microburst  detection  and  gust 
front  detection,  are  in  general  worldng  well  Minor 
software  changes  are  being  made  to  inprove  system 
performance. 


Issues  surrounding  building  accurate  ground 
clutter  suppression  maps  became  apparent  and  will 
need  to  be  addressed  at  future  sites.  Additional 
testing  is  scheduled  for  Spring  1993  after  which  final 
results  will  be  published  in  the  5th  International 
Aviation  Weather  Conference  to  be  held  in  August, 
1993. 
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QUANTIFYING  AIRPORT  TERMINAL  AREA  WEATHER  SURVEILLANCE 

REQUIREMENTS* 


Marilyn  M.  Woifson  and  Cynthia  A.  Meuse 
Lincoln  Laboratory,  Massachusetts  Institute  of  Technology 
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1.  INTRODUCTION 

The  Federal  Aviation  Adminisuation  (FAA)  Terminal 
Area  Surveillance  System  (TASS)  research,  engineering,  and 
development  program  was  initiated  in  part  to  address  future 
weather  sensing  needs  in  tire  terminal  area.  By  the  early  21st 
century,  planned  systems  such  as  theTerminal  Doppler  Weath¬ 
er  Radar  (TDWR)  and  Airport  Surveillance  Radar-9  (ASR-9) 
will  be  well  into  their  designed  life  cycles.  Any  new  terminal 
weather  surveillance  system  should  be  designed  to  address  ex¬ 
isting  deficiencies.  Key  unmet  weather  sensing  needs  include 
detections  of.  Uaie  3-dimensional  winds  (vs.  radial  compo¬ 
nent),  winds  m  the  absence  of  precipitation,  wahe  vortices,  to¬ 
tal  lightning,  bail,  icing  conditions,  clear  an  turbulence,  haz¬ 
ardous  weatliercells(wilh  adequatetime  andspacerescli'tion), 
cloud  cover  and  cloud  bases  (including  layers),  fog,  and  visi¬ 
bility  (Runway  Visual  Range),  as  well  as  predictions  of.  the  at¬ 
mospheric  conditions  mentioned  above,  wind  shifts,  micro- 
bursts,  tornadoes,  and  snow/rainfall  rates  (Evans  1991a, 
McCarthy  1991). 

In  this  paper,  we  investigate  the  premise  that  hazardous 
weather  cells  arc  not  currently  being  measured  with  adequate 
time  and  space  resolution  m  the  terminal  area.  Since  anew  sur¬ 
veillance  system  should  be  based  on  knowledge  of  storm  dy¬ 
namics,  we  have  performed  a  preliminary  study  of  update  rate 
(using  rapid  scan  radar  data  collected  in  Orlando),  and  spatial 
resolution  required  to  detect  rapidly  developing  tliunder- 
stonns  and  precursors  to  the  low  altitude  hazards  such  as  rm- 
crobursts  that  theyproduce.  Other  aspectsof  afutuiciadarsys- 
tem  such  as  multi-parameter  techniques  required  to 
discriminate  betwceniccandwaterphaseprecipitation, etc.  are 
not  considered. 

2.  APPLICATIONS  FOR  TERMINAL  AREA 

CONVECTIVE  WEATHER  INFORMATION 

Past  studies  have  shown  tliatweatheris  the  primary  cause 
of  serious  air  traffic  delay  at  the  nation’s  m.ajor  airports  (Weber 
et  al.  1991)  and  a  recent  study  at  Lincoln  Laboratory  has 
shown  that  this  delay  may  even  be  underestimated  wiUi  the 
current  reporting  system.  Thunderstorms  and  heavy  fog  ac¬ 
count  for  tlic  largest  fraction  of  weatlier  related  delay.  Some  of 
tliis  delay  can  be  nutigatedby  more  accurate  detection  and  pre- 
dictionofweatherimpactcdfliglitroutes.allowingefficientre- 
routing  to  take  place.  Tlr  is  will  provide  an  important  econormc 
and  safety  benef it  in  tlie  future,  as  increased  air  traffic  demands 
maximization  of  capacity  at  our  existing  airports. 

The  newly  deployed  TDWRs  .ind  ASR-9s  will  provide 
a  significant  increase  in  safety  m  tlie  terminal  area,  but  these 
systems  were  not  specifically  designed  to  reduce  weather  re¬ 
lated  aviation  system  delays.  Tlie  FAA  has  recently  mitiatcu 
theIntcgratcdTcmiinalWealherSystem(rrWS)program,bfc- 
mgdcvclopedatLincolnLiboratory,to  pros  idc  improved  avi¬ 
ation  weather  informatioii  m  tlie  termuial  area  by  integratmg 

•Thtt  A’Ofk  waijponsofMbyibe  r«5<r»l  A>iU'onAdmmiJlraiion.TK»view»c/prM>«J 
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craineni. 


data  and  products  from  tiie  FAA  .ind  NWS  sensors  (Evans 
1 99  lb).  A  key  objective  of  tlie  ITWS  program  is  to  increase  tlie 
effective  airport  acceptance  rate  in  adverse  weatlier  by  provid¬ 
ing  information  to  support  tlie  Terminal  Air  Traffic  Control 
Automation  (TATCA)  program  TATCA  tools  increase  tlie  ef¬ 
ficiency  of  individual  controller  tasks  and  provide  a  dynamic, 
overall  plan  for  traffic  management  tJiroughout  die  terminal 
control  region  (Andrews  and  Welch  1989).  Tlius,  reliable  .anal¬ 
yses  and  forecasts  of  weatlier  impacted  air  routes,  runway 
availability,  and  clear-air  winds  for  direr,  t  use  by  TATCA,  traf¬ 
fic  managers,  and  pilots  arc  a  major  goal  of  tlie  ffWS  system 
Tlie  design  of , my  future  terminal  weather  surveillance  system 
will  have  to  consider  tliis  close  coupling  of  weather  sensing 
andforecasting,andairtrafficc<ipacityandefficicnc>  "nhance- 
ment  programs. 

The  first  deployment  of  tlie  ITWS  system  will  inch’de  a 
Microburst  Prediction  product,  and  may  also  include  a  storm 
location  prediction  based  on  projection  of  existing  storms  ac¬ 
cording  to  correlation  tracking  information .  Tlicsc  predictions 
are  performed  on  existingdctectable  radar  reflectivity  regions 
and  thus  are  necessarily  short  term.  Longer  term  predictions  of 
weather  impacted  airspace  will  require  predictions  of  storms 
that  ha  ve  yet  to  develop,  and  tlius  will  citlier  depend  heavily  on 
heuristic  niles  for  convective  initiation  (Mueller  and  Wilson 
1993)  or,  perhaps  more  likely  because  of  accuracy  require¬ 
ments,  will  require  gridded  4-D  data  assimilation  -numerical 
forecasting  techniques  such  as  those  being  developed  at 
NOAA’s  Forecast  Systems  Laboratory  (Slicrretz  1991).  Tliis 
type  of  gridded  analysis  system  is  already  being  used  in  3-D 
for  an  ITWS  terminal  winds  representation  (Wilson  et  al. 
1993).  Witli  or  without  gridded  numerical  forccasimg  tech¬ 
niques,  convective  forecasts  will  be  requ  ired  in  a  fu  11 60  km  ra¬ 
dius  region  around  the  airport  at  least. 

3.  UPDATE  RATE 

Of  all  the  unmet  terminal  weadier  sensing  needs,  the  de¬ 
sire  to  predict  rapidlydevelopingconvective  weather  is  a  driv¬ 
ing  factor  in  the  proposed  required  update  rale  for  a  TASS 
weather  sensor  Keeler  (1991)  assigns  a  critical  time  scale  of 
0  5  -  1 .0  min  for  thunderrtorms,  .ind  1.0-  2.0  min  for  wide¬ 
spread  rain.  Tliere  is  evidence  that  a  1  min  update  rate  is  re¬ 
quired  for  or  til  anderstorm  outflow  detection.  Study  has  dem¬ 
onstrated  the  need  for  a  similar  update  rate  for  tliundcrstomi 
life  cycle  predictions  (Carbone  et  al.  1985). 

To  investigate  tliis,  rapid  scan  iiicasurcmcnts  were  made 
witli  tlie  IDWTI  testbed  radar  operated  by  Lincoln  L.iboratory 
in  Orlando,  FL,  where  the  typically  very  unstable  env  uonment 
leads  to  rapidly  developing  and  decaying  tliundcrstorms.  A 
special  TASS  scan  was  designed  to  cover  a  complete  volume 
in  1  min  witli  a  set  of  elev.ition  .ingles  conipar.iblc  to  a  tme 
TDWR  hazardous-weatlier  mode  scan  (sec  Fig,  1).  Tlie  scan 
provides  uniform  coverage  up  to  1 4  km  AGL  for  a  storm  at  20 
km  range.  Suiuibly  long  rapid  scan  cLitasctson  multiple  micro¬ 
burst-producing  storms  were  gatliercd  on  Aug.  5  ,  Sept.  22, 
and  Sept.  26 1992.  For  tlie  identified  cells  (5  total).  p.iramcters 
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Figure  1.  The  TASS  Rapid  Volume  Scan  consisis  of  10 
elevation  angles  (0  5, 3  7. 7.3, 10  9. 14.7. 18  6, 22.5, 26  6. 
30  7, 35  0°  ).is60°  wide,  and  takes  1  min  to  execute  with  the 
TDWR  testbed  radar. 


such  as  theaverage  vertically  mtegrated  liquid  waier(\TL)and 
ihe  heiglil  of  the  center  of  mass  (CM)  were  computed.  Trends 
in  these  parameters  arc  used  in  the  ITWS  .Vlicrobursl  Predic¬ 
tion  Algorithm  (Wolfson  et  al.  1 993)  to  identify  growing  thun¬ 
derstorms  and  to  predict  their  collapsing  phase,  which  leads  to 
microburst  wind  shear  at  tlic  surface.  These  parameters  mea- 
suredevery  nunute  arc  compared  with  tlie  identical  parameters 
derived  from  hybrid  TDWR  volume  scans  made  up  of  3  TASS 
scans  each,  with  the  lowest  elevation  angles  coming  from  the 
first  TASS  scan,  tlic  middle  angles  from  the  second,  and  the 
highest  angles  from  tlic  third.  As  an  example,  tlic  VILdata  for 
Aug  5  are  shown  in  Fig  2a.Tlie3-min  VIL  lags  and  peaks  lat¬ 
er  tlian  the  1  -min  VIL,  which  shows  much  more  detailed  fluc- 
tuations.The  height  of  center  of  mass  data  show  a  similar  pat¬ 
tern  (Fig.  2b). 


TIME  (GMT) 

Figure  2.  (a)Plot  of  VII,  computed  for  a  storm  cell  scanned 
with  the  TASS  Rapid  Volume  Scan  on  August 5, 1992,  and  of 
VIL  com puted from  3-niinh\  hr  id  scans  made from  the  same 
data.  The  surface  outflow  AV  derived from  the  1 -min  scans 
IS  also  shown.  The  black  circle  alAV=15  indicates  the  mi- 
croburst  onset  and  the  white  circle  at  A  V=I0.  wtndshear  on¬ 
set.  (b)  As  in  (a)  but  for  height  of  center  of  mass  (CM). 


Weusedthe  prototype  ITWS  Microburst  Predictionalgo- 
ritlim  to  quantify  tlie  advantage  of  tlie  TASS  scan  ( 1  min)  over 
tlie  hybrid  TD'W  scan  (3  mm).  In  the  algoritlim.  a  region  is 
first  idcp"'"ed  and  tracked  based  on  a  significant  nse  in  VIL. 
among  oilier  features.  Tins  rise  must  persist  for  2  volumes  for 
a  track  to  be  established.  If  tlie  region  has  been  tracked  twice, 
a  prediction  can  be  issued  as  soon  as  a  drop  in  center  of  mass 
is  detected,  assuming  tlie  quantity  of  VBL  present  at  tliat  time 
is  sufficient  to  produce  a  niicroburst-slrength  outflow.  Table  1 
shows  tlie  achievable  prediction  lead  times  for  tlic  TASS  and 
hybrid  TDWR  scans.  Tlie  TASS  strategy  allows  an  .average  -3 
mm  gie.iter  prediction  lead  iinic.  extciidmg  tlie  TDWR  aver¬ 
age  of  2-4  min  for  tliese  (weak)  cases  to  5-6  min  (Tlicse  re¬ 
sults  may  change  as  tlie  Microburst  Prediction  algontlim 
evolves  from  its  prototypical  to  its  final  form ) 

4.  SPATUL  RESOLUTION 

Higli  spatial  resolution  is  required  to  detect  microbursts, 
tornadoes,  etc.  at  very  low  elevation  angles,  but  it  may  be  pos¬ 
sible  to  trade  resolution  for  a  more  rapid  update  rate  at  upper 
levels.  To  investigate  tins,  we  created  a  model  stomi  ellipsoid 
10  km  higli.  6  km  wide,  witli  a  central  core  reflectivity  of  65 
dBZ  at  5  km  AGL,  winch  decreased  linearly  to  tlie  outer  edges 
of  tlie  ellipsoid.  We  compared  tlic  TDWR  scan  to  an  exper¬ 
imental  low  resolution  fan  beam  sc.an  consisting  of  6  beams, 
each  5  ^  in  elevation,  centered  at  5, 1 5, 25, 35, 45,  and  55°.  By 
ushg  only  6  broad  beams  to  scan  tlic  volume  instead  of  1 1  nar¬ 
row  beams  (TDWR),  tlieupdatcratecouldtheoretically  be  im¬ 
proved. 

Since  VILand  CM  measurements  arc  crucial  tothe  ITWS 
Microbursi  Prediction  algorithm.  iJtcse  parameters  as  a  func¬ 
tion  of  range  were  compared  Tlie  ellipse  was  moved  in  range 
from  0  to  25  kni,  and  ■‘scaimcd'’evcry  km  with  botli  strategics. 
Figure  3  shows  tliat,  altliougli  tlic  cone  of  silence  over  the  ra¬ 
dars  leads  to  incorrect  values  inside  of  4  km  range,  tlic  values 
foi  VO  strategies  arc  not  dramatically  different  at  near 
r.r  •’thcTDWRVILfallssbonoftlictruc  VILbecause 

intcij  on  cannot  recreate  unsampled  peak  reflectivity  re¬ 
gions.  ‘d  the  1  km  influence  r.idius  c.-ui  include  distant  low 
VIL  values.)  At  longer  range,  our  studies  have  shown  that  tlie 
large  variability  of  the  fan  beam  VIL  is  due  to  the  gaps  in  eleva¬ 
tion.  and  tlie  large  rise  in  fan  beam  CM  is  due  tothe  rising  center 
height  of  tlic  1 5°  beam,  which  dominates  the  CM  calculation. 
For  microburst  prediction,  the  apparent  changes  in  CM  witli 
rangehavefargrcatcrpotcntialforcausingfalsealamistliando 
tlie  changes  in  VIL. 

Since  future  algoritlims  and  especially  any  numerical 
forecasting  techniques  will  undoubtedly  use  gridded  r.idar 
data,  we  also  investigated  the  effect  of  Cartesian  grid  resolu- 

Table  1.  Microburst  prediction  lead  times  for  5  cells 
scanned  with  TASS  scan  and  hybrid  TDWR  scan.  For  cases 
with  2  entries,  first  is  lead  time  for  onset  of  10  mis  outflow, 
second  for  15  mis. 


Date 

TDWR  lead 

TASS  lead 

Outflow  DV 

1992 

(min) 

(min) 

(m/s) 

August  5 

3-6 

6-9 

17.5 

September  22A 

5 

7 

11.6 

September  22B 

3 

5 

14.8 

September  22C 

0  -  1 

2-3 

19.3 

September  26 

0 

6 

Average  10  m/s: 

2.2 

5.2 

5  cases 

Average  15  m/s: 

3.5 

6.0 

2  cases 
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Distance  from  Radar  (km) 

Figure  3.  Apparent  VILand  height  of  CM  of  model  ellipsoid 
scanned  with  TDWR  and  fan  beam  strategies. 


lion.  We  discovered  that  there  was  very  little  difference  be¬ 
tween  0. 1  km  and  0.5  km  for  tlie  TDWR  scan  (Fig.  4).  At  1 .0 
km.  the  measured  parameters  were  still  very  close  to  theu-  0. 1 
km  values,  but  the  vanabihty  of  CM  with  range  increased.  Tlic 
ITWS  prototype  Microburst  Prediction  algoritlini  currently 
uses  1  km  resolutionfaltliougliO  Skmvertical  resolution  is  be¬ 
ing  considered),  and  the  gridded  Terminal  Winds  product  uses 
2  km  horizontal  resolution  at  50  mb  heiglit  intervals. 

5.  SUM.VIARV  AND  CONCLUSIONS 

We  have  considered  here  the  key  unmet  needs  in  Uie  air¬ 
port  terminal  area  of  adequate  coverage  both  m  tinicandspacc 
for  detection  and  prediction  of  hazardous  weather  cells.  Be¬ 
cause  tliundcrstorms  evolve  on  such  rapid  time  scales,  tins  re¬ 
quirement  may  be  a  driver  for  tlic  design  of  any  new  termuial 
weatlicr  surveillance  system  It  is  our  contention  that  require- 


Figure  4.  Apparent  height  of  center  of  mass  of  model  ellip¬ 
soid  scanned  with  TDWR  strategy,  and  gridded  at  0.1, 0.5. 
and  1 .0  hn  resolution.  A  I  0  km  interpolation  radius  was 
used  at  each  resolution. 


ments  for  update  rate  and  spatial  resolution  must  be  developed 
in  the  context  of  likely  use  of  tlie  data  -  in  endded  prediction 
systems,  data  assimilatingnumcrical  models,  and  for  specific 
improvements  in  fliglit  route  plaiuiuig.  wcatlier  avoidance, 
and  terminal  air  traffic  capacity  ."uid  efficiency. 

We  have  developed  a  metliodology  towards  quantifying 
the  wcatliersurveillancc  requirements  forupdate  rate  .uid  cov¬ 
erage  of  convective  storms  By  makingspecial  radar  measurc- 
menis  at  higli  volume  update  rates,  .and  comparing  tlicm  witli 
derived  lower  update  rate  scans,  die  benefit  witli  respect  to 
somedetcctioncrprcdictionalgontlimormodcl  can  be  quanti¬ 
fied.  Likewise,  by  sampling  idealized  or  actual  high  resolution 
numerical  model  data  of  storms,  trade-offs  on  required  cover¬ 
age  resolution  can  be  made.  Based  on  Uiese  studies  and  oilier 
research  performed  at  Lincoln  Laboratory  for  die  TDWR  and 
ITWS  programs,  we  can  draw  several  conclusions. 

1.  Wliole  volume  coveragetoat  least  4kni  altitude  in  the  ter¬ 
minal  area  is  required  every  1. 0-2.5  min.  The  actual  selected 
update  rate  will  depend  on  the  perceived  cosifbcncfit  ratio 
St.atc-of-die-art  algoridims  and  numerical  models  should  be 
used  to  calculate  diis  ratio 

2.  A  “cone  of  silence",  i.e.  an  unscanned  ''cgion  over  die  ra¬ 
dar.  is  not  acceptable  over  die  airport. 

3.  Surface  updates  of  0.5-1 .0  min  are  required  over  the  air¬ 
port  and  approach/departure  corridors  for  timely  outflow 
detection. 

4.  Higli  resolution  data  (100-250  m  range  gate  spaemg. 
1-1.5“  azimmlial  resolution)  arc  required  m  die  boundary  lay¬ 
er,  i  e.  die  lowest  -2  kni.  Contiguous  beam  coverage  is  desir¬ 
able  since  windspeedchangesrapidlywidiheiglitncarilie  sur¬ 
face.  It  IS  cmcial  to  dcterrmiiingnew  areas  of  convection  dial 
boundary  layer  forcing  (convergence)  be  adequately  mea¬ 
sured.  and  diis  requires  low  altitude  coverage  even  at  the  far 
ranges  of  the  tcmiinal  area. 

5.  Above  -2  kni.  resolution  requirements  can  be  relaxed. 
Wider  beam  widdis  can  be  tolerated,  but  the  desired  resolution 
at  long  range  must  still  be  taken  into  account. 
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A  MACHINE  INTELLIGENT  GUST  FRONT  ALGORITHM  FOR  DOPPLER 

WEATHER  RADARS 


Richard  L.  Oclanoy  and  Seth  W.  TVoxcl 

Lincoln  Laboratory,  Massachusetts  Institute  of  Technology 
Lexington,  MA  02173-9108  USA 


1.  INTRODUCTION 

Gust  fronts  generated  by  tliunderstorms  can  seriously  affect 
the  safety  and  efficiency  of  airport  operations.  Lincoln  Lab¬ 
oratory,  under  contract  wiih  the  Federal  Aviation  Adminis¬ 
tration  (FAA),  has  bad  a  significant  role  in  the  development 
of  two  Doppler  radar  systems  tliat  are  capable  of  detecting 
low-altitude  wind  shears,  mcluding  gust  fronts,  in  the  air¬ 
port  tenniiul  control  area.  These  systems  are  the  the  latest 
generation  Airport  Surveillance  Radar,  cnltanced  with  a 
Wind  Shear  Processor  (ASR-9  WSP)  and  the  Terminal 
Doppler  Weather  Radar  (TDWR). 

Gust  fronts  produce  signatures  that  arc  observable  to  varying 
degrees  in  reflectivity  and  Doppler  velocity  data  generated 
by  these  radars.  In  Doppler  velocity  images,  gust  fronts  arc 
recognizable  as  zones  of  converging  winds.  In  reflectivity 
images,  gust  fronts  appear  as  thin  lines  of  increased  intensity, 
which  occur  as  the  result  of  rain,  dust,  insects,  or  debris  be¬ 
ing  lofted  and  concentrated  at  die  leading  edge  of  the  front. 
An  existing  automated  gust  front  detection  and  forecasting 
algoritlira,  developed  principally  for  TDWR  data  and  called 
in  this  paper  the  Advanced  Gust  Front  Algorithm  (AGFA), 
has  achieved  respectable  levels  of  performance  using  these 
data  (Eilts,  1991  and  Merritt,  1989).  With  clear,  unambigu¬ 
ous  radar  signatures  AGFA  performs  reasonably  well.  The 
challenge  is  in  constructing  an  algorillim  tliat  c.an  detect  mar¬ 
ginally  detectable,  ambiguous  cases  witliout  incurring  unac¬ 
ceptable  false  alarm  rates. 

Several  sources  of  ambiguity  e.xist.  For  example,  gust  front 
thin-line  signatures  can  be  obscured  by  large  areas  of  preci¬ 
pitation.  Velocity  convergence  signatures  can  vanish  when 
gust  front  orientations  result  in  bad  Doppler  viewing  angles. 
Gust  front  signatures  can  also  be  mimicked  by  other  natural 
phenomena,  such  as  flocks  of  birds,  clouds  of  dust  stirred  up 
at  construction  sites,  elongated  low-intensity  rain  echoes, 
and  ground  clutter.  Finally,  gust  fronts  can  have  very  low 
radar  cross-section  densities,  sometimes  below  the  sensitiv¬ 
ity  of  the  radar  system,  making  detection  difficult. 

The  ASR-9  WSP  provides  a  less  expensive  alternative  to 
the  TDWR  as  a  terminal  weather  radar  (Weber,  1989).  Al¬ 
though  not  originally  intended  for  weather  imaging,  this 
fai>-bcam  Doppler  radar  generates  images  of  sufficient  qual¬ 
ity  that  gust  fronts  can  be  identified  and  tracked.  However, 
versions  of  AGFA  ad-ipted  for  ASR-9  WSP  data  have  per¬ 
formed  poorly.  The  primary  reason  for  the  lack  of  perfor¬ 
mance  is  the  reduced  gain  and  lowered  sensitivity  inherent 

•The  work  described  here  was  sponsored  by  the  Federal  Avi¬ 
ation  Adnunistntion.  The  United  States  Government  assumes 
no  liability  for  its  content  or  use  thereof. 


in  tlie  fan-beam  design  of  the  ASR-9.  With  lowered  sensi¬ 
tivity.  clear  air  velocity  estimates  are  unreliable.  Conse¬ 
quently.  few  wind  convergence  signatures  arc  visible,  forc¬ 
ing  AGFA  to  rely  on  only  the  thin-line  portion  of  its 
algorithm.  Unfortunately,  the  reduced  sensitivity  also  makes 
faint  thm-line  signatures  more  fragmented  and  harder  to  re¬ 
solve  from  background. 

The  radars  themselves  are  sufficient  for  tlx:  task  of  gust  front 
detection,  since  experienced  human  observers  can  detect  arxl 
track  gust  fronts  in  images  generated  by  these  radars.  And 
yet,  sufficiently  high  detection  rates  with  few  false  alarms 
has  been  an  elusive  goal  for  developers  of  automated  gust 
front  detection  algorithms.  Skilled  human  intcTX’cicrs  rely 
upon  spatial  and  temporal  contextual  information  and  assim¬ 
ilate  weak,  uncertain,  ambiguous,  and  even  contradictory  ev¬ 
idence.  Humans  arc  also  adept  at  conditionally  fusing  in¬ 
formation  from  various  sources,  reflecting  knowledge  that 
different  signatures  can  have  varying  reliability  that  depends 
on  situational  context.  In  contrast,  such  traits  of  perceptual 
intelligence  have  been  notoriously  and  surprisingly  difficult 
to  implement  in  computer  vision  systems. 

2.  ALGORITHM  DESCRIPTION 

A  general-purpose  approach  to  object  recognition,  which 
was  originally  developed  in  the  context  of  automatic  target 
recognition  (ATR),  has  been  incorporated  in  a  Machine  In¬ 
telligent  Gust  Front  Algoritlim  (MIGFA).  Use  of  tire  term 
"machine  intelligence”  in  particular  reflects  tlic  use  of  two 
new  techniques  of  knowledge-based  signal  processing. 

Tire  first  technique,  functional  template  correlation  (FTQ 
(Dclanoy.  1992)  is  a  generalized  matched  filter  incorporating 
aspects  of  fuzzy  set  tlicory.  For  comparison,  standard  2-D 
cross  correlation  relies  upon  a  kernel  tliat  is  essentially  a  sub- 
image  consisting  of  expected  image  values.  In  contrast,  the 
kernel  of  a  functional  template  consists  of  a  set  of  integers 
tliat  each  correspond  to  a  unique  scoring  function.  Each 
scoring  function,  given  an  image  value  as  input,  returns  a 
score  reflecting  how  well  that  image  value  matched  expecta¬ 
tions  for  a  given  location  on  the  kernel.  The  results  of  all 
scoring  functions  within  the  functional  template  are  then  av¬ 
eraged,  resulting  in  a  score  in  tlie  range  [0,1]. 

By  increasing  or  decreasing  the  interval  over  which  affirm¬ 
ing  scores  (i.e.,  >  0.5)  arc  relumed,  scoring  functions  can 
encode  varying  degrees  of  uncertainty  with  regard  to  what 
image  values  arc  allowable.  But  in  addition,  knowledge  of 
how  a  feature  or  object  appears  in  sensor  imagery  can  be  en¬ 
coded  in  scoring  functions.  With  various  design  strategies, 
the  intcrfeiing  effects  of  occlusion,  distortion,  noise,  and 
clutter  can  be  minimized.  As  a  consequence,  matched  fillers 
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customized  fo''  specific  applications  using  FTC  are  generally 
more  robust  tlian  traditional  signal  processing  operations. 
Hie  output  of  FTC  is  a  map  of  values  in  tbc  range  [0,1], 
each  of  which  reflects  a  belief  that  tile  shape  or  object  im¬ 
plicitly  encoded  in  a  functional  template  is  present  at  that 
image  location.  In  our  ATR  systems,  FTC  has  been  used 
primarily  as  a  direct,  one-step  means  of  3-D  object  detec¬ 
tion  and  extraction.  In  MIGFA,  FTC  is  used  for  edge  detec¬ 
tion,  tliin-line  filtering,  tliin-liiie  smoothing,  shape  match¬ 
ing,  and  lioniotopic  tliinning  of  shapes. 

The  second  major  tool  is  tlie  >'se  of  "interest"  as  a  medium 
for  data  fusion  and  for  assimilating  evidence  at  the  pixel  lev¬ 
el  (Dclanoy,  1991).  An  interest  image  is  a  map  of  evidence 
(values  in  the  range  {0,1])  for  some  feature  that  is  selectively 
indicative  of  an  object  being  souglii  (note  that  the  output  of 
FTC  can  be  an  interest  image).  Higlicr  pixel  values  reflect 
greater  confidence  tliat  the  intended  feature  is  present  at  that 
location  Given  the  assumption  tliat  the  output  of  any  fea¬ 
ture  detector  can  be  configured  as  an  interest  image,  evi¬ 
dence  from  any  number  of  registered  sources  of  information 
can  be  easily  combined  using  simple  or  arbitrarily  complex 
rules  of  aritlinKtic  or  fuzzy  logic.  Clusters  of  high  values  in 
tlic  combined  interest  image  arc  tlicn  used  to  guide  selective 
attention  and  serve  as  tlie  input  for  object  extraction.  In 
practice,  we  often  use  several  weakly  or  inconsistently  dis 
criminiiting  feature  detectors  to  mutually  support  or  compen¬ 
sate  for  each  otlicr,  resulting  in  relatively  robust  perfor¬ 
mance. 

The  sy.stcm  block  diagram  in  Fig.  1  illustrates  tlie  configura¬ 
tion  of  tlic  ASR-9  WSP  version  of  tlie  system.  Input  images 
V  (Doppler  velocity  image)  and  (tZ  (reflectivity  image)  arc 
passed  to  six  simple,  uidcpendent  feature  dctceton  that  use 
FTC. 


omcTows 


Figure  1.  MIGFA  block  diagram. 


deviation)'.  A  second  feature  dctccto,"  labelled  SEX-MO¬ 
TION,  apphes  a  thin  line  filter  to  the  difference  of  two  se¬ 
quential  reflectivity  images.  Since  differencing  suppresses 
those  signals  that  are  stationary,  tlic  SD-MOTION  detector 
tends  to  highlight  only  those  thin-lines  tliat  are  moving. 
The  third  feature  detector,  labelled  OUT-OF-TRIP,  uses 
FTC  to  identify  range  ambiguous  echoes  of  naore  distant 
weather.  Areas  believed  to  reflect  out-of-trip  wcatlier  arc 
given  low  interest  values.  The  next  two  feature  detectors, 
TD-DZ  and  DZ-MOTIOM  perform  tlie  same  kind  of  thin- 
line  analysis  ss  is  done  for  TL-SD  and  SEX-MOTION.  Tlie 
final  feature  detector.  ANTICIPATION,  will  be  discussed 
below.  The  outputs  of  tlic  several  feature  detectors  arc  com¬ 
bined  as  a  weiglilcd  average  to  form  a  combined  interest 
image. 

From  tlie  combined  interest  image,  fronts  arc  extracted  as 
cliains  of  points.  The  chains  extracted  from  a  radar  scan,  col¬ 
lectively  called  an  event,  are  integrated  wi'Ji  the  prior  history 
by  establishing  point-to-point  correspondence.  Heuristics 
arc  applied  at  tliis  point  to  reject  chain  points  that  have  ap¬ 
parent  motion  tliat  is  improbable.  Tbc  updated  history  is 
used  to  make  predictions  of  where  points  along  the  front  will 
be  located  at  some  future  time.  Such  predictions  arc  used  in 
tlic  processing  of  subsequent  images,  specifically  in  tlie  fea¬ 
ture  detector  called  ANTIQPATION,  which  places  liigli  in¬ 
terest  values  wherever  fronts  arc  expected  to  be  and  by  so 
doing  selectively  sensitizing  tlic  j^stem  to  detecting  gust 
fronts  at  specific  locations.  Anticipation  is  tuned  fo  tliat  it 
will  not  by  itself  automatically  trigger  a  detection,  but  when 
averaged  w  ith  oilier  uitcrcst  images  will  support  weak  evi¬ 
dence  tl--  would  otherwise  be  insufficient  to  trigger  a  detec¬ 
tion. 

A  second  version  of  the  gust  front  algorillim  has  been 
consUuctcd  for  TDWR  data,  differing  only  in  the  set  of  fea¬ 
ture  detectors  used.  Differences  in  tlie  TDWR  detector  set 
prini.Trily  reflect  tlic  greater  dependence  on  Doppler  data  for 
Tuiding  .areas  of  convergence. 

3.  RESULTS 

A  test  set  of  ASR-9  WSP  data  collected  in  Orlando,  Florida 
during  AGFA  field  testing  in  1991,  contains  9  different  gust 
fronts  in  a  set  of  450  images  (IS  hours).  A  human  interpret¬ 
er  looking  at  tlic  same  ASR-9  WSP  data  detected  280 
instances  of  the  9  gust  fronts  tracked  by  tlic  radar.  Four  fig¬ 
ures  of  merit  arc  shown  for  each  of  the  two  algoritlims.  Hie 
probability  of  detection  is  tlic  number  of  detections  made  by 
each  algorithm  as  a  percent  of  bunian-dclcctcd  instances  of 
gust  fronls.  The  probability  of  false  alarm  (PFA)  is  the  num¬ 
ber  of  false  alarms  divided  by  tlic  total  number  of  algorithm 
de'.cctions.  In  addition  to  simply  identifymg  fronts,  the  hu¬ 
man  interpreter  delimited  the  length  of  each  detected  front. 
Detection  quality  was  further  assessed  by  comparing  tlie 
length  of  tlic  front  as  estimated  by  each  algorithm  against 

1.  Because  measurements  within  gust  front  thin  line  have 
higher  signal  to  noise  ratios  than  the  surrounding  clear  air  mea¬ 
surements,  gust  fronts  show  up  in  SD  as  thin  lines  of  lower  sig¬ 
nal  variance. 


The  first  one,  marked  Tl^SD,  uses  an  FTC  filter  for  thin- 
lincs  in  a  map  of  local  EXopplcr  variance,  called  SD  (standard 
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that  indicated  by  the  human  interpreter.  Tlie  percent  Icngtlr 
detected  ”'LD)  is  the  length  detected  expressed  as  a  percent 
of  the  ^tb  delimited  by  the  human  interpreter.  The  percent 
of  fili/>  cflgtb  detected  (PFD)  reflects  situations  where  Uie 
detecterl  gust  front  lengilis  extended  beyond  what  the  human 
interprete'  could  see. 

Table  1  compares  performance  of  MIGFA  against  tlic  pre¬ 
viously  constructed  AGFA,  which  uses  more  conventional 
methods  of  signal  processing  and  computer  vision.  Tlie  first 
two  columns  indicate  tliat  MIGFA,  relative  to  AGFA,  sub¬ 
stantially  increased  the  POD  while  decreasing  the  PFA. 
Similarly,  the  PLD  reflects  the  improvement  in  detection 
rate.  However,  the  increased  PFD  (from  13%  to  33%) 
would  .suggest  that  MIGFA  was  doing  a  worse  job  of  dis¬ 
criminating  the  extent  of  individual  fronts.  In  order  to  better 
understand  why  MIGFA  was  extending  fronts  beyond  what 
the  human  interpreter  believed  appropriate,  we  rescored 
AGFA  and  MIGFA  results  against  human  interpretations  of 
TDWR  data  taken  at  the  same  time  as  tlie  ASR-9  WSP  data. 
Wj  assume  that  truth  derived  from  the  more  reliable  TDWR 
data  is  more  accurate  than  that  for  ASR-9  data.  Tlicse  re¬ 
sults,  shown  in  Table  2,  confirm  the  general  trend  of  the  first 
3  figures  of  merit  shown  in  Table  1 .  However,  now  the  PFD 
is  quite  low,  essentially  tlie  same  as  Jiat  for  AGFA.  Conse¬ 
quently,  these  results,  along  wiili  analyses  of  individual 
cases,  leads  us  to  believe  tliat  the  MIGFA  delected  fronts 
were  in  fact  more  accurate  tluan  detections  made  by  tlie  hu¬ 
man  interpreter  given  tlie  same  ASR-9  WSP  data. 


Table  1.  AGFA  and  MIGFA  performance  on  ASR-9 
WSP  data  as  scored  against  human  interpretations. 


Gust  Fronts 

Gust  Front  Length 

POD 

PFA 

PLD 

PFD 

AGFA 

56.7 

46 

38.9 

12.9 

MIGFA 

88.1 

0.6 

86.2 

33.4 

Table  2.  AGFA  and  MIGFA  performance  on  ASR-9 
WSP  data  as  scored  against  human  interpretations  of 
matching  TDWR  data. 


Gust  Fronts 

Gust  Front  Length 

POD 

PFA 

PLD 

PFD 

AGFA 

42.6 

3.2 

21.0 

4.2 

MIGFA 

75.1 

0.0 

58.7 

6.4 

This  same  versi.''n  of  MIGFA  was  installed  on  an  ASR-9 
WSP  in  Orlando,  Florida  hr  operational  testing  during  tlie 
summer  of  1992,  tlie  results  of  which  arc  shown  m  Table  3. 
During  the  period  from  1  August  to  20  September,  MIGFA 
correctly  detected  and  tracked  approximately  75%  of  all  gust 
fronts  identified  by  human  intr  'prcicrs  examining  ASR-9 
WSP  data.  TDWR  and  anemoir.ctcr  data  were  also  used  for 
verification.  Those  gust  fronts  that  were  missed  either  bad 


reflectivity  levels  near  the  sensitivity  Lmils  of  the  ASR-9  or 
were  obscured  by  storm  cells  along  the  front.  The  false 
alarm  rate  was  under  2%.  Although  TDWR-based  truth  has 
net  yet  been  compiled  for  tlicse  data,  an  analysis  of  individu¬ 
al  cases  again  indicates  that  tlie  relatively  higli  PFD  (21%) 
consists  largely  of  believable  extensions  of  gust  fronts  tliat 
were  not  identified  by  tlie  human  interpreter. 


Table  3.  Results  of  MIGFA  operational  testing  on 
ASR-9  WSP  data  collected  in  Orlando  during  sum¬ 
mer  1992. 


Gust  Fronts 

Gust  Front  Lcngtli 

POD 

PFA 

PLD 

PFD 

MIGFa 

75.4 

1.8 

80.8 

21.1 

4.  STATUS 

MIGFA  represents  a  substantial  improvement  in  perfor¬ 
mance  over  pr.  'ions  efforts  and  is  llic  prime  candidate  for 
deployment  in  production  ASR-9  WSP  systems.  Tlie 
TDWR  version  of  MIGFA  is  scheduled  for  testing  in  sum¬ 
mer  1993.  Adapt.ations  of  tlie  techniques  used  in  MIGFA 
arc  currently  being  used,  or  arc  being  considered,  for  otlicr 
weatlicr  detection  problems,  including  microburst  prediction 
and  sensor  fusion. 
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1.  INTRODUCTION 

Tlie  Federal  Aviation  Administration  (FAA)has  been  spon- 
soruig  Lincoln  Laboratory  in  its  effort  to  develop  and  test 
wcatlicrdetectionalgoritlimsforthe Terminal  Doppler  Weath¬ 
er  Radar  (TDWR).  An  automated  microbursl  detection  algo- 
ritlim(Merrittetal.,  1989)operateson  the  TDWRradial  veloc¬ 
ity  data  and,  based  on  the  shear  and  velocity  difTcFcncc  along 
tlic  radial,  outputs  regions  which  arc  hazards  to  aviation.  This 
algorithm  has  been  operating  since  1987  in  Denver,  Kansas 
City,  and  Orlando  and  is  part  of  tlic  operational  TDWR  being 
deployed  across  tlie  country.  One  issue  which  continues  to 
cause  concent  for  automated  windshcar  detection  is  micro¬ 
burst  asynunctry.  Asymmetry,  or  aspect  angle  dependence,  in 
microbursls  refers  to  outflows  which  have  a  divergent  surface 
outflow  slrengtli  or  extent  tliat  varies  depending  on  the  view¬ 
ing  angle  of  the  radar. 

The  TDWR  is  a  single-Doppler  radar,  tlierefore,  an  asym¬ 
metric  microburst  may  be  underestimated  or  go  undetected  if 
tlie  radar  is  viewing  the  event  from  an  aspect  angle  where  the 
strengtli  of  tlie  outflow  is  weak.  Past  work  by  Wilson  et  al. 
(1984),  Eilts  (1987,1988),  and  Hallowell  (1990)has  indicated 
that  some  microbursts  arc  liiglily  asymmetric.  Strengtli  asym¬ 
metries  (maximutn/minimuni  strengtli  over  all  viewing 
angles)  from  these  past  studies  ranged  from  1 .3  to  as  higli  as 
6  0.  Hallowell  (1990)  using  Denver  data  examined  27  Denver 
microbursts  (96observations)  and  found  strength  asymmetries 
from  1.3  to  3.8  with  a  median  of  1.9.  However,  this  previous 
work  has  been  limited  in  scope  to  Denver  and  Oklaiioma 
(plains)  microbursts,  and  may  have  used  assumptions  about 
tlie  data  which  introduce  false  or  apparent  asymmetry. 

2.  APPARENT  ASYMMETRY 

Previous  investigators  selected  dual-Doppler  microburst 
events  using  the  following  criteria: 

•  Intersection  angle  of  beams  between  30°  and  150° 

•  Tilt  limes  of  the  two  radars  within  one  minute 

•  Elevation  angles  of  botli  radars  less  than  1.0° 

hi  tlic  courscof  studying  microbursl  events,  we  have  found 
lliat  while  tlicse  assumptions  are  valid  for  large  scale,  slow 
moving  and  developing  wind  fields,  llicy  are  not  sufficiently 
strict  formicrobuist  outflows.  Weutilized  simulated  three-di¬ 
mensional  velocity  data  of  a  symmetric  microbursl  obtained 
usingthc  WME  (Wisconsin  Model  Engine)  sub-cloud  model 
(Anderson,  1992).  Radial  data  were  extracted  from  lire  simu¬ 
lated  data  from  twohypotlietical  radars  situated  15  kmfroralhc 
microburst  center  at  various  elevations  angles,  beam  widths, 
limes,  and  intersection  angles.  We  llicn  input  lliese  lilts  in  vari¬ 
ous  combinations  to  gatiicr  data  on  how  Uie  apparent  asymme¬ 
try  of  tlic  event  changed  for  each  parameter.  Using  these  model 


•The  work  described  here  was  spon^r cd  by  the  Federal  AMatioa  Adminisua* 
don.  TYte  United  States  Government  assumes  no  liability  for  lU  content  or  use 
(hereof. 


results  we  were  able  to  quantify  tlie  effects  of  various  dual- 
Doppler  coordination  parameters  on  tlic  measured  asymmetry 
of  tlie  real  microburst  data.  Tlie  followmg  base  configuration 
was  used  (changes  arc  noted  in  each  sub-section).  Fictitious 
radars  located  15  kmfrom  microburst  center  at  90°  beam  inter¬ 
section  angles,  1.0°  beam  width,  0.4°  elevation  angles,  and  a 
gate  spacing  of  150  meters. 

2.1.  Diial-Donoler  Process  Itself 

The  wind  synthesis  process  which  operates  on  the 
smootlied  polar  radial  velocity  data  and  creates  a  two-dimen¬ 
sional  Cartesian  windfield  grid  intioduces  asymmetry  of  its 
own.  By  extracting  two  tills  of  fictitious  radar  data  at  90° 
angles  and  keeping  all  oilier  parameters  tlic  same,  we  find  a 
dual-Doppler  microburst  field  which  yields  .an  asymmetry  of 
1.04.  This  apparent  asymmetry  has  not  been  removed  in  the 
graphs  that  follow. 

2.2.  Horizontal  Beam  Intersection  Aiiples 

The  angle  of  intersection  between  tlic  radar  radial  data  isex- 
tieniely  important  in  determining  tlic  quality  of  the  dual-Dop¬ 
pler  wind  field.  By  holding  all  oUicr  parameters  consl.int  and 
changing  only  the  intersection  angle  of  tlie  radars,  we  see  from 
Fipre  1  that  intersection  angles  less  llian  45°  yield  increasing 
apparent  asymmetry  results.  _ _ 
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Figure  1.  Apparent  Asymmetry  vs.  hitcrscction  /\nglc 

2.3.  Temporal  Vari.itions 

Tlie  time  between  tilts  is  .'mother  key  measure  of  data  quali¬ 
ty,  in  microbursts  one  minute  can  mean  a  change  in  strengtli  of 
over  10  ni/s.  Looking  at  tlie  observations,  we  fmd  that  90%  of 
tlie  Iminutechangesindiffcremial  velocity  arc  lessilian7  m/s. 
The  microbursl  model  used  has  a  peak  differential  velocity  of 
51  mis  and  one  minute  prior  to  tlic  peak  a  velocity  of  44  m/s. 
By  matching  extracted  tills  from  different  15  sec  time  steps 
(witli  all  other  parameters  tlie  same),  we  find  tliat  data 
compared  with  45  secs  or  less  lime  differential  yields  very  low 
apparent  asymmetries  (Figure  2).  Care  should  be  taken  in  ap¬ 
plying  tills  to  all  microbursts;  rapidly  developing  or  decaying 
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Figure  2.  Apparent  Asymmeu-y  vs.  Time  Difference 
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microbursts  will  be  more  of  a  problem  even  in  a  30  sec  time  in¬ 
terval. 

2  4.  Ream  Wiclili.  F.levaiioii  Ancle  Differences 

On  examinmg  tlicse  parameters  separately  it  is  difficult  to 
discern  a  pattern.  Wlial  we  fmally  found  was  that  the  extent  of 
vertical  beam  overlap  was  tlie  mipuitant  feature,  not  a  specific 
beam  width  or  elevation  angle.  By  processing  one  tilt  at  0.0® 
elevation  anditstiltpairatuicreasmglyhigherelevationangles 
we  fmd  beam  ovcrlapto  be  tlieovemdiiigconcem  for  apparent 
asynunotiy  fFigurc  3).  Tlie  simulauon  used  for  tliis  analysis 


_ Beam  Overlap  (%) _ 

Figure  3.  Apparent  Asymmetry  vs.  Beam  Overlap 


_ Fraction  of  Maximum  Strengtli _ 

Figure  4.  Vertical  Strengtli  Profile:  Simulated  (solid)  & 
Measured  (dashed') _ 

was  for  a  very  strong  microburst,  and  its  vertical  profile  is 
steeper  than  tliat  found  by  Biron  and  Isamingcr  (1991)  for 
Denver  nticrobursts.Figure4showstlie  vertical  profijeofboth 
tlic  simulated  data  (solid  line)  and  the  Denver  observations 
(dashed  line)  Wiiile  tlie  overall  simulated  profile  presents  a 
worst  case  scenario,  the  lowest  200  m  are  fairly  comparable, 
and  tins  is  where  tlie  30%  andgreateroverlapanalysis  wasper- 
fomied. 

2.5.  Recommendations 

Tlie  followmgcriteriaweicuscdin  tliis  analysis, and  should 
be  used  by  otlier  researchers,  to  limit  the  affects  of  apparent 
asynmictry  on  asymmetry  analyses: 

•  Intersection  angles  between  45°  and  135° 

•  Tilt  lime  differences  <  30  secs,  and 

•  Radar  beam  overlap  by  50%  or  more,  or  beam  centers 
wiili  in  50  m  (ill  is  guarantees  lliat  each  radar  catches  at 
least  a  portion  of  ilic  otlier  radar’s  center  beam). 

3.  DATA 

Tlicdalau.ed  were  collected  in  Denver,  CO  (1988),Kansas 
City,  MO  (1989)  and  Orlando,  FL  (1990).  At  each  site,  there 
were  two  radars  operating'  tlie  MIT/LL  TDWR  testbed  radar 
(FL-2)  and  tlie  University  of  North  Dakota  C-Band  Doppler 
Radar  (UND)  FL-2  and  UND  boUi  have  one  degree  beams. 


FL-2  was  an  S-Band  radar  prior  to  1990  but  modified  to  C- 
Band  in  1990  Tlie  radar  scanning  was  coordinated  to  cover 
niicrobursts  which  occurred  in  favorable  dual-Doppler  re¬ 
gions.  For  each  scan  of  an  event,  the  two-dimensional  surface 
wmd  field  was  calculated  using  tlie  multiple  Doppler  radar 
synthesis  system  suggested  by  Brown  el  al.  (1981).  The  raw 
two-dimensional  wmd  fields  were  then  smootliedonceusiiig 
a  simple  3-by-3  median  filler,  wiili  4  of  9  points  required  to 
be  valid. 

Microbursis  were  selected  by  examining  tlie  two-dimen¬ 
sional  windficld  for  divergence  regions.  This  subjective  ex¬ 
amination  piocess  yielded  over  1000  observations  and  some 
100  events  (multiple  observations  of  tlie  same  microbursl). 
Some  observations,  taken  at  the  begmning  or  end  of  an  event, 
were  not  true  outflows  and  were  removed  from  tlie  analysis  by 
dropping  observations  which  had  minimum  differential  velo¬ 
cities  less  Ilian  3.5  m/s.  Wiili  this  restriction,  a  total  of  859  ob¬ 
servations  of  microbursts  were  examined  for  asymmetry;  a 
breakdown  by  site  is  shown  in  Table  1. 

Table  J. Breakdown  by  Site  of  Data  Examined 


Days 

Events 

Observations 

Denver 

6 

36 

476 

Kansas  City 

7 

27 

163 

Orlando 

4 

22 

220 

Totals 

17 

85 

859 

4.  METHODOLOGY 

Tlie  meiliodology  used  to  analyze  microburst  events  from 
dual-Doppler  is  essentially  the  same  as  that  detailed  by  Hallo- 
well  (1990).  Briefly,  an  analyst  examines  the  wind  field  and 
draws  a  polygon  around  a  microburst  outflow.  Gridpoints  of 
velocity  data  witliin  tlie  polygon  are  ilien  intercompared  tocal- 
culate  velocity  differences.  Every  gridpoint  pair  has  its  own 
orientation  angle  wiili  respect  to  Norili  (based  on  a  fictitious  ra¬ 
dar  15  kni  from  ilic  event,  see  Figure  5),  from  which  weean  de- 
tcmiinc  what  the  peak  differential  is  at  various  orientation 
angles 

5.  RESULTS 

Three  sites  and  859  observations  of  microbursis  were  ex¬ 
amined  for  tliis  study.  Wlicn  tlie  strict  criteria  from  Section  2.5 
arc  applied,  we  Fmd  a  dcrmiic  shift  toward  lower  asymmetry  as 
shown  in  Figure  6.  The  median  asymmetry  ratio  for  the  rc- 
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slrictedsci  ofdaia  (498  observations)  is  1 .78  whilcall  oilier  ob¬ 
servations  have  a  median  asymmetry  over  2.0  (361  observa¬ 
tions).  By  no  means  has  all  llic  apparent  asymmetry  been 
eliminated;  resiriciingcasesiorcmove  ail  apparentasymmetry 
would  leavenoobservalions  to  work  wiili.  However,  a  reason¬ 
able  estimate  of  the  lowest  possible  real  asymmetry  curve  can 
be  made  based  on  die  known  apparent  asymmetry  error  left  in 
the  data.  The  measured  asymmetry  can  be  obtained  by  apply¬ 
ing  ilie  following  formula; 

^measured”  Afcai*A(,nie  (liff*A(xa(n_iinele*Aovcrtap*Adual 

This  formula  assumes  a  worst  case  scenario,  where  the 
asymmetries  due  to  temporal  variations  (Aume  dill),  beam  in¬ 
tersection  angles  (Abeam.angie).  elevation  angle  overlap 
(Aovcriap).  and  ilic  dual-Doppler  process  itself  (Ajuai)  com¬ 
pound  one  another  Sometimes  die  apparent  asymmetries  may 
be  oriented  such  dial  diey  actually  counteract  one  anodier,  al¬ 
though  the  test  data  studied  indicates  diis  is  less  likely.  By  re¬ 
stricting  die  data,  I  have  attempted  to  limit  the  impact  of  each 
apparent  asymmetry  to  under  5%  for  A<jual.  Adme.diff. 
Abcam.aneie.  and  10  under  15%  for  Aovertap  By  dividing  the  re¬ 
stricted  measured  asymmetry  curve  by  our  estimated  error 
(1.15*1  05^-1  33),Iobiainameasurcofdice.xpecu:dcumula- 
tive  distribution  of  real  asymmetry  (Figure  7)  Tlie  chart  indi¬ 
cates  the  real  asymmetry  median  could  be  as  low  as  1.34,  but 
the  actual  answer  likely  lies  between  the  two  curves. 

If  we  examine  asyiunieiry  on  asitc  by  sue  basis,  wc  find  that 
foi  each  site  die  restricted  data  set  yields  lower  asymmetries 
dian  unrestricted  data  (not  shown).  Figure  8  shows  the  cumula¬ 
tive  distribution  of  asymmetry  ratios  for  each  of  die  diree  sites 
examined.  Orlando  and  Denver  turn  out  to  have  very  similar 
distributions  of  .asymir  clry,  widi  medians  of  1 .72  and  1 .76, rc- 
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Figure  7.  Cumiiblivc  Frequency  of  Ratios  for  Restricted  D.iM. 
Measured  (grey)  and  No  Apparent  Asymmetry  (black). 


speclively.  Kansas  City  shows  a  tendency  low.ardhiglier  asym¬ 
metry  (median  1 .9),however  1989  was  an  atypical  climatolog¬ 
ical  year  diere.  Because  ofdiis,  die  quality  and  quantity  of  the 
data  collected  in  Kansas  City  may  have  been  msufficieni  to 


make 

any  firm  conclusions  about  midwest  niicrobursis. 
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Figure  8.  Cumukilive  Frequency  of  Ratio.s  for  Restricted  Data: 

Denver  (grey),  Kansas  Cily  (thin  black),  and  Orlando  (thick  black). 

6.  CONCLUSIONS 

Previ3usstudicsofmicrobursiasymmciryliavenotconsid- 
eied  die  impact  of  radar  configurations,  ilicreby,  introducing 
apparent  asymmetry  mto  die  .analysis.  By  usingsimulateddata 
we  are  able  to  estimate  die  impact  of  temporal  difference,  verti¬ 
cal  beam  overlap,  and  horizontal  beam  mterseclion  angle  fac¬ 
tors  in  processing  and  analyzing  dual-Doppler  microburst 
events.  A  total  of  859  microburst  observations  were  exanuned 
from  three  geographical  regions.  VVe  find  dial  overall  asymme¬ 
try  distribudons  arc  lower  dian  haa  been  found  m  all  previous 
studies,  and  that  differences  m  asynimeiry  between  sites  such 
as  Orlando  and  Denver  are  muiinial  Overall,  the  measured 
asymnieuy  ratios  m  die  observation  data  vary  from  1 . 1  to  4.0 
(with  less  dian  5%  over  3.0)  and  have  a  nicdi.Tn  value  of  1  78. 
In  addition,  estimatmg  and  removing  die  residual  errors  of  ap¬ 
parent  asymmetry  from  die  observation  data  set  yields  a  dis¬ 
tribution  of  I.O  to  3.0  and  a  median  of  1.34. 
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1.  INTRODUCTION 

Microbursts  are  small  scale,  low  altitude  wind 
shea/  phenomena  which  have  been  associated  with  many 
recent  aircraft  accidents.  Microbursts  arise  ftomthunder- 
stormsandarecharacterizedbyintense  downdraftswhich 
spread  out  after  impacting  the  earth’s  surftice  and  display 
strong  divergent  outflows  of  wind.  Th^  are  usualty 
associated  with  heavy  rainfall  (Wolfton,  1988). 

The  Tferminal  Doppler  WeatherRadar(TWDR) 
program  was  the  first  attempt  at  microburst  detection  with 
a  ground-based  radar  in  the  airport  terminal  area.  Im¬ 
proving  safety  was  its  primary  goal,  and  test  operations  in 
Denver,  Kansas  Gty,  and  Orlando  have  shown  it  hi^ly 
successful  in  identifyingmicrobursts.  In  general,  this  iden¬ 
tification  hasbeen  performed  with  a  >  90%ProbabDityof 
Detection  (POD)  and  a  <  10%  Probability  of  False  Alarm 
(PFA)  (Campbell  and  Olson,  1987). 

The  Integrated  Tferminal  Weather  System 
(iTWS)  seeks  to  enhance  this  ability.  Microburst  Predic¬ 
tions  will  be  produced  in  addition  to  detections  (Wolfton 
et  al.,  1993).  A  miaoburst  trend  product,  giving  predic- 
tionsofincreasingraicroburstintensityalongrunway  corri¬ 
dors  over  short  time  periods  (2-3  minutes),  is  to  be 
introduced.  Thismicroburst  trend productwill  involve  the 
ability  to  predict  the  future  location,  size,  and  intensity  of 
the  microburst  It  is  largely  due  to  the  microburst  trend 
product  thatthephilosophybehindmicrobuistdetectionis 
being  revised. 

Although  the  TDWR  algorithm  is  successful  in  detecting 
the  area(s)  of  hazard,  the  output  representation  is  not 
suited  for  tracking  microbursts.  The  ITWS  algorithms  at¬ 
tempts  to  alleviate  this  by  providing  one  output  shape  for 
eachdowndraft.  TheTDWRalertingisfundamentallyloss 
based,  that  is,  the  severity  of  the  hazard  is  indicated  by  the 
strengthofthssurfacedivergencccouplet.  However,ifthis 
divergence  is  not  over  a  small  area,  an  aircraft  will  experi¬ 
ence  little  or  no  performance  deficit  The  ITWS  algorithm 
captures  this  information  by  examining  the  divergence 
shear  (rate  of  change  in  velocity)  as  well  as  the  loss. 
Ground-based  dopplerradar  observation  hasshown ,  and 
instrumented  aircraft  penetrations  have  confirmed  (Mat¬ 
thews  and  Berke,  1993,  Campbell  et  al.,  1992).  that  the 
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shearwithinamicroburstishighlynonuniform.  ThelTWS 
algorithm  will  use  an  additional,  imbedded  warning  shape 
to  indicate  especially  hazardous  regions  of  a  microburst 

This  paper  erqrlains  the  initial  designof  the  rrWS 
microburst  detection  algoiithm  and  illustrates  some  early 
results.  Thefinalsectkmconcentratesonthetestingplans 
for  algorithm  testing  and  the  intended  enhancements  to  its 
capabilities. 

Z  ALGORITHM  DESIGN 

The  algorithm  is  comprised  of  four  primary  ele¬ 
ments:  radial  shear  calculation,  segment  formation,  region 
formation,  and  alarm  generadotL  The  algorithm  is  intended 
to  find  microbursts,  and  weaker  wind  shear  events,  out  to  a 
30.0  km  range,  but  3S.0  km  is  actually  processed  to  alleviate 
edge  effects.  Surface  scans  from  the  TDWR  are  available 
qjproximately  cveryminute.  TherrWSmicroburstalgorithm 
v^l  process  dl  scans  equal  to,  or  below  1.0  degrees  elevation 
angle. 

The  radial  shear  is  calculated  as  the  spatial  deriva¬ 
tive  of  the  radial  velocity  field.  The  base  velocity  data  is  first 
median  filtered,  using  a  range  adaptive  filter  size  (approxi¬ 
mately  1km  by  1km),  and  a  least  squares  fit  of  a  line  segment 
to  the  data  is  performed  as  in  Figure  1.  The  slopeofthefitline 
segment  is  written  as  the  radial  shear  at  that  gate.  Both  the  me¬ 
dian  filtering  and  the  regression  fitting  output  a  valid  value 
only  if  at  least  a  certain  percentage  (currently  50%)  of  base  ve¬ 
locity  values  are  valid. 


Figure  1.  Illustration  of  shear  calculation.  The  slope  of  the 
heavy  line  is  assigned  as  the  shear  value. 

The  segment  forming  module  processes  the  radial 
shear  map  to  locate  contiguous  segments  along  a  radial  above 
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«setofus«-defined  thresholds.  Currently,  two  thresholds,  5.0 
m/s/km  (low  shear)  and  10.0  m/s/km  (high  shear)  are  used. 
The  segment  must  meet  a  minimum  length  to  be  considered 
valid.  Each  valid  segment  is  then  extended  until  either  (a)  the 
average  shear  along  its  length  falls  below  the  threshold,  (b)  a 
gate  withanegative  shear  value  is  encountered,  or  (c)  an  exces¬ 
sive  number  of  consecutive  invalid  shear  values  are  found.  In 
this  way,  a  segment  is  not  considered  unless  a  minimum  shear 
is  present,  but  the  length  of  the  segment  is  not  rigidly  tied  to  the 
points  above  thre*  old.  The  loss  (velocity  difference)  across 
the  segment  is  calculated  over  a  different  scale,  found  as  the 
maximum  segment  length  while  the  average  shear  remains 
above  2.5  m/s/km,  or  until  conditions  (b)  or  (c)  are  found. 


Range  — ► 


Figure  2.  Depiction  of  Segment  forming  process.  The  final 
segment  endpoint  I  was  placed  by  condition  (a),  and  endpoim 
n  was  placed  by  condition  (b). 

The  Region  formation  module  processes  the  set  of 
shear  segments  at  each  threshold  level  into  regions  of  high 
shear.  Regions  are  built  by  associating  segments  on  adjacent 
radials  which  overlap  by  a  user  defined  percentage.  A  circul  ar 
shape  is  fit  to  the  segments  of  the  region  using  an  optimization 
technique.  The  parameters  of  the  circle  (x-center,  y-center, 
radius)  are  altered  to  minimize  the  sum  of  the  minimum  dis¬ 
tance  from  each  segment  end  to  the  circle  perimeter  and  the  dis¬ 
tance  of  the  segment  end  to  the  circle  center.  The  optimization 
is  unidimensional  (each  parameter  is  optimized  one  at  a  tii.'te), 
since  it  was  determined  that  the  extra  computation  for  a  multi¬ 
dimensional  technique  was  not  necessary.  The  region  is  dis¬ 
carded  if  a  minimum  number  of  segments  is  not  iiKluded  or  if 
the  area  of  the  best  fit  circleis  below  a  threshold.  Regions  from 
the  high  and  low  shear  threshold  levels  are  associated  together 
if  the  percent  area  of  intersection  is  higher  than  a  user-defined 
amount.  Those  high  shear  regions  without  an  accompanying 
low  shear  region  are  discarded. 


The  low  shear  regions  are  passed  through  a  final 
alarm  generation  test.  Ifthemaximum  loss  from  any  segment 
in  aregion  is  greater  than  30  knots  and  the  peak  shear  within 
the  region  is  greater  than  10.0  m/s/km,  a  microburst  alert  is 
generated.  If  themaximum  loss  exceeds  15  knots  and  thepeak 
shear  is greater  than5.0m/s/km,awindshear  alert  isgenerated. 
All  high  shear  regions  associated  with  the  alarm  regions  are 
also  output  for  display. 

3.  RESULTS 

The  algorithm  has  thus  far  been  executed  on  12  mi¬ 
croburst  cases,  using  TDWR  testbed  data  collected  from  Or¬ 
lando  (8  cases),  Kansas  City  (2  cases),  and  Denver  (2  cases). 
On  this  data  set  it  has  demons  tratedthatitcanatleastmatchthe 
performance  of  the  TDWR  algorithm  in  ideruifying  micro- 
bursts  and  quantifying  their  intensi^.  Addidoital  studies  are 
necessary  to  ensure  that  the  false  alarm  rate  is  also  comparable. 

The  additional  goals  of  a  more  accurate  downdraft 
identification  and  output  shape  consistency  as  an  aide  to  track¬ 
ing  have  produced  mixed  results .  Isolatrxi  events,  such  as  that 
shown  in  Figure  3, 4,  and  5  are  successfully  characterized  as 
originating  from  single  downdrafts.  As  aresult,  the  algorithm 
is  successful  in  assigning  event  labels  and  following  a  consis- 
tenttrackoverthelifetimeoftheevenL  Situations  with  multi¬ 
ple,  in  teractingdowndrafts  are  curren  tlyidentifiedmoreerrati- 
cally.  This  originates  from  an  improper  identification  of  the 
downdrafts,  andresultsindifiiculttointerpretapparerumerges 
and  splits  of  events.  Recent  evidence  about  the  correlation  be¬ 
tween  peak  shear  locatioirs  and  downdrafts,  and  additional 
aloft  information  which  will  be  integrated  from  thelTWS  Mi¬ 
croburst  Prediction  algorithm  (Wolfson  et.  al.,  1993)  indicate 
that  this  problem  can  be  alleviated  in  the  short  term. 


Figure  3.  The  loss  and  peak  shear  values  over  the  course  of 
an  Orlando,  1990  miaobur.st,  as  determined  by  the  ITWS  Mi¬ 
croburst  Detection  algorithm. 
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Figure  4.  The  location  track  of  the  same  microburst  event  as 
in  Figure  3.  The  center  of  the  circle  ou^t  by  the  ITWS  algo¬ 
rithm  is  used  as  the  location,  and  the  truth  has  been  determined 
by  hand  by  an  expert  meteorologisL 

4.  FUTURE  DIRECTION 

The  IfWS  algorithm  will  be  tested  in  real-time  in 
the  ITWS  testbed  site  in  OrlarKio,  Florida  from  5/93  through 
9/93 .  This  is  expected  to  give  us  a  large  dataset  for  the  inevita¬ 
ble  modification  and  Cne  tuning  of  the  algorithm.  The  ITWS 
Microburst  Prediction  and  Microburst  Trend  algorithms  will 
also  be  running  in  Orlando,  begirming  7/93 .  This  will  help  ex¬ 
amine  the  strong  interactions  between  the  algorithms. 

There  are  several  enhancements  to  the  algorithm 
which  are  expected  to  be  made  in  thenext  calendar  year.  Other 
sensors  ate  to  be  examined,  particularly  the  Low  Level  Wnd 
Shear  Alert  System  (LLWAS),  for  integration  with  the  ITWS 
algorithm.  This  will  ensure  that  a  consistent  single  alert  is 


made  for  the  enure  termmidl  area.  Attempts  will  be  made  to 
compensate  for  microburst  asymmetry  and  adjust  for  aircraft 
altitude  in  microbuist  intensity  warnings.  Aloft  irtformation, 
of  the  kind  the  ITWS  Microburst  Prediction  algorithm  is  con¬ 
sidering  (Wolfson  et.  al.,  1993),  will  be  used  to  provide  a  more 
accurate  idendficadon  of  the  strength  and  locadon  of  the 
downdraft.  This  is  important,  since  the  verdcal  winds  are  not 
capnrred  by  low  elevadon  angle  ground  based  radar  scans,  and 
the  downdraft  is  also  a  source  of  significant  hazard  for  an  air¬ 
craft. 
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Figure  5.  Comparison  of  (a)  TDWR  and  (b)  ITWS  algorithm  shapes  for  time  20;  17  in  the  event  of  Figure  3  and  4.  A  high  shear 
region  is  cross-hatched  within  the  outer  ITWS  shapte. 
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1.  INTRODUCTION 

In  the  past  decade,  a  great  deal  of  effort  has  been  in¬ 
vested  in  developing  ground-based  wind  shear  detection  sys¬ 
tems  for  major  U.S.  airports.  However,  there  has  been  a  lack  of 
research  in  devclopingaquantitativerclationshipbetween  the 
wind  shear  iiazards  detected  by  ground-based  systems  and  the 
actual  hazard  experienced  by  an  aircraft  flying  through  the  af¬ 
fected  air  space.  To  date,  tlie  mam  thrust  of  tlie  verification  ef¬ 
forts  for  ground-basedsystemshasbeen  to  ensure  that  thesys- 
tem  accurately  detect  and  report  the  presence  of  the 
meteorological  phenomena  that  cause  potentially  hazardous 
windshear.  There  is  a  subtle,  but  potentially  important  differ¬ 
ence  between  detecting  tlie  presence  of  a  microburst  and  de¬ 
tecting  the  presence  of  an  aviation  hazard.  Wlh  this  in  mind, 
it  would  seem  prudent  to  rigorously  determine  what  correla¬ 
tion  exists  between  the  wind  shear  warnings  that  are  generated 
from  ground  systems  and  tlie  performance  impact  on  aircraft 
flying  througli  the  impacted  airspace.  The  operational  demon¬ 
stration  of  the  tested  Terminal  Doppler  Weather  Radar 
(TDWR)  in  Orlando,  Florida  along  witli  thetestingof  airborne 
Doppler  radar  systems  created  a  unique  opportunity  to 
compare  extensively  tlie  ground-based  windshear  reports 
with  in-situ  aircraft  raeasurcmenis. 

Tliis  paper  presents  the  results  from  69  microburst 
penetrations  flown  in  1990  and  1991  by  the  University  of 
Nortli  Dakota (UND), tlie  National  AcronauticsandSpace  Ad¬ 
ministration  (NASA)  Langle  y  Research  Center,  and  Rockwell 
Collins  under  surveillance  of  the  Lincoln-operated  TDWR 
testbed  radar  The  primary  goal  of  the  research  was  to  deter¬ 
mine  tlie  relative  accuracy  of  several  methods  designed  to  gen¬ 
erate  a  numerical  microburst  hazard  index,  called  the  F  factor, 
from  ground-based  Doppler  radar  data.  It  is  hoped  that  this 
work  will  provide  both  a  qualitative  and  quantitative  basis  for 
the  discussion  and  assessment  of  microburst  hazard  reporting 
for  ground-based  microburst  detection  systems. 

Die  Integrated  Airborne  Wind  Shear  Program  is  a 
joint  NASA/Federa!  Aviation  Administration  (FAA)  program 
with  the  objective  to  provide  tlie  technology  base  tliat  will  per¬ 
mit  low  altitude  wuidshear  risk  reduction  through  airborne 
detection ,  warn  ing,  and  avoidance.  Additionally,  the  program 
aims  todcmonslratc  tlie  practicality  and  ulilityof  real-time  as¬ 
similation  and  synthesis  of  ground-derived  windshear  data  to 
supportexecutivelevelcockpitwaniingandcrcw-centeredin- 
formation  display.  Lincoln  Laboratory  joined  this  effort  and 
provided  tlie  weatlicr  radar  ground  support  and  some  of  the 
*The  work  described  here  was  sponsored  by  the  Na¬ 
tional  Aeronautics  and  Space  Administration.  The 
United  States  Government  assumes  no  liability  for  its 
content  or  use  thereof. 


post-fliglit  data  analysis  for  NASA's  microburst  penetration 
flights  in  Orlando,  Florida. 

2.  F  FACTOR  EQUATION 

Currently,  ground-b.ased  and  airborne  windshear 
sensors  characterize  tlie  microburst  hazard  in  two  different 
ways.  Airborne  systems  characterize  the  hazard  in  terms  of  F 
factor(Bowlcs,  1^0),  a  volumetric  parameter  Uiat  measures 
the  rate  of  change  of  aircraft  energy  Ground-based  systems 
such  as  TDWR  report  the  strengtli  of  llic  windshear  event  as  a 
single  number  representing  peak  point-to-point  loss.  There¬ 
fore,  it  was  necessary  to  process  the  radar  bfise  data  to  eompute 
the  F  faetor  for  each  penetration.  This  was  accomplished  by 
two  methods.  Die  first  metliod  involved  using  the  TDWR  al- 
goritlim  loss  estimate  and  the  second  involved  a  shear  calcula¬ 
tion  from  the  radial  velocity  data.  Both  methods  used  a  Lui- 
cota  modified  version  of  the  F  factor  equation  proposed  by 
Roland  Bowles  of  NASA  Langley  Rc.ccarcii  Center  (Bowles. 
1990): 

"r  =  +  i)  =  +  F-  « 

where  K’  is  a  constant,  A  V  is  the  velocity  difference,  AR  is  the 
distance  overwhich  tlie  vclocitydiffercncc  occurred.  GS  is  tlie 
groundspreed  of  the  aircraft,  and  h  is  the  hciglit  of  Uie  radar 
beam.  The  F  factor  equation  is  composed  of  two  terms,  tlie 
horizontal  term(Fh,  effect  ofheadwind/tail  wind  loss  on  tliean- 
eraft),  and  tlie  vertical  tcrm(Fv,  effect  of  tlie  downdraft  on  air¬ 
craft  performance).  The  Doppler  radar  is  only  capable  of  mea¬ 
suring  tlie  wind  component  along  a  radial,  tlicrefore,  the 
aircraft  generally  flew  fliglit  paths  along  a  radial.  Diis  com¬ 
pensated  for  inconsistences  between  radar  and  aircraft  mea¬ 
surements  in  the  horizontal  term,  however,  some  metliod  was 
needed  to  estimate  tlie  vertical  term  from  tlie  Doppler  data. 
This  was  done  by  employing  a  simplified  mo  i  of  tlie  mass 
continuity  equation.  Die  outflow  region  is  vu  as  a  cylm- 
der  with  the  hciglit  of  the  radar  beam  acting  as  tlie  lop  of  tlie 
cylinder.  Because  tlie  ground  acts  as  acap  preventmg  flow  out 
the  bottom,  what  flows  into  die  top  (i.e.  downdraft)  must  exit 
through  tlie  sides  of  the  cylmder  as  horizontal  outflow.  Witli 
this  r^odel,  the  outflow  is  directly  proportional  to  the  down¬ 
draft  uependmg  upon  the  radius  of  the  cylmder  and  the  decel¬ 
eration  profile  of  tlie  downdraft  witli  heiglit. 

3.  ESDMATTON  OF  F  FACTOR  USING  TDWR 
ALGORDDM  OUTPUT 

The  TDWR  microburst  algorithm  defines  a  micro- 
burst  outflow  region  by  fitting  a  race  track  shaped  icon  around 
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Figure  1  TDWR  algorithm  output  vs 
aircraft  total  F  factor 


groups  of  radar  radial  velocity  segments  tliat  show  sufficient 
shear  along  tlteir  length.  Site  adaptable  parameters  control  the 
maMmum  sl^e  of  a  shape,  the  nimimum  shape  radius,  and  otli- 
cr  shape  ch  ..ractcrislics.  Each  icon  is  assigned  a  value  denotuig 
the  pcak-to-pcak  velocity  difference  contained  within  tlie 
shape.  This  ranged  hur « die  largest  to  the  second  largest  seg¬ 
ment  AV  dcpendtngupon  tlie  number  of  s.’gmcnts  contained 
within  the  shape.  This  oeak-to-peak  velocity  difference  was 
taker,  to  be  tlie  AV  term  of  equation  (1).  and  the  AR  term  was 
based  upon  tlie  85th  percentile  shear  within  the  shape. 

For  each  of  Uie  69  microburst  penetrations,  an  F  fac¬ 
tor  was  calculated  from  the  output  of  tlie  Lincoln  version  of  the 
TDWR  nucroburst  algonthm  using  the  techniques  described 
above.  Figure  1  is  a  plot  of  tlie  TDWR  estimated  total  F  factor 
(TDWR  Ft),  compared  to  tlie  in  situ  F  factor.  From  the  figure 
it  can  be  seen  that  the  computed  TDWR  Fp  was  consistently 
higher  than  tlie  tit  situ  F  factor.  Most  notable  is  that  the  estima¬ 
tion  was  biased  especially  high  for  the  NASA  events. 

Tlie  unexpectedly  high  values  of  TDWR  Ft  may  be 
due  to  several  factors,  the  most  obvious  of  which  is  tliat  the  F 
factor  computed  from  the  rmcroburst  alarms  assumes  tliiit  the 
shear  li.as  a  constant  value  at  all  points  witliin  the  alarm’s 
boundary.  Tliis  assumption  is  incorrect,  and  since  the  aircraft 
sampled  on  ly  a  smal  I  portion  of  tlie  area  enclosed  by  the  micro- 
burst  alarms,  it  is  quite  possible  that  on  many  occasions  they 
missed  the  localized  "hotspot”  of  shear  that  caused  the  large 
value  reported  by  tlie  TD\W.  Tlic  test  pilots  indicated  during 
interviews  tliat  tliey  occasionally  avoided  tlie  most  severe  por¬ 
tion  of  a  storm  intentionally  due  to  fliglit  safety  considerations. 


4.  ESTIMATION  OF  F  FACTOR  USING  TDWR 
SHEAR  MAP 

F  factors  canbe  computed  from  the  TDWR  testbed’s 
radial  velocity  data  by  creating  a  nwp  of  tlie  radi<il  shear.  To 
do  til  is,  tlie  rad.v  base  data  were  first  subjected  to  a  data  quality 
editing  process  and  then  velocity  dc.ali.asing.  The  data  quality 
editing  consisted  of  clutter  removal,  point-target  editing,  and 
range  obscurauon  editing  Next,  tlie  velocity  field  was  median 
filtered  using  a  slidmg  wmdow  of  approximately  500  meters 
X  500  meters.  The  actual  radial  shear  computation  for  each 
range  gate  was  m.ade  by  performing  a  Ic.ist  squ-ircs  fit  on  seven 
gates  centered  about  tlie  point.  Witli  die  TDWR  radar’s  150 
meter  gate  spacmg,  diis  resulted  m  a  fit  over  a  radial  distance 
of  1050  meters.  This  general  metliod  of  shear  computation  is 


Aircraft  F  Factor 
Figure  2  TDWR  shear  map  vs  air¬ 
craft  tola!  F  factor 


Similar  to  work  done  by  Britl(  1 992)  for  NASA’s  airborne  Dop¬ 
pler  windshear  detection  system. 

The  F  factor  can  tJien  be  estimated  along  die  trajec¬ 
tory  of  the  aircraft  by  using  die  c  losest  radial  shear  value  as  cal¬ 
culated  from  the  TDWR  base  data  Figure  2  shows  the  peak 
total  F  factor  as  estimated  from  the  TDWR  shear  map  versus 
the  peak  total  in  situ  F  factor  A  comp.arison  widi  Figure  1 
shows  improved  agreement  between  die  TDWR  and  in  situ  F 
factor  .at  dieexpense  of  an  increased  incidence  of  underestima¬ 
tion.  Examination  of  the  significant  citses  of  underestimation 
reveals  that  in  nearly  all  cases  die  error  was  due  to  die  aircraft 
encountering  a  large  downdraft  diat  was  not  predicted  by  die 
shear  map  F  factor  calculation. 


5  COMPARISON  OF  AIRCRAFT  AND  RADAR 
ESTIMATES  OF  HORIZONTAL  F  FACTOR 


As  mentioned  in  Section  3,  it  is  useful  to  look  at  die 
horizontal  and  vertical  terms  of  die  F  Factor  when  atiemptmg 
to  analyze  the  success  and  failure  of  the  various  csiim.aiion 
tecliniqucs.  From  equation  1 ,  the  horizontal  temi  can  becalcu- 
lateddnectly  from  die  TDWR  shear  mapdata  us  jig  die  follow¬ 
ing  formula  (Note.  K’  is  equal  to  one  because  d.e  shear  m.ap  is 
a  one  kilometer  shear): 


Fh 


=  AV/GS\ 
AR\  g  / 


(2) 


Figure  3  compares  thchorizontal  termofdicF  fac¬ 
tor  as  estimated  from  die  TDWR  shear  map  and  die  aircraft, 
nic  shear  mapprovides  a  fairly  good  estinwte  of  die  horizon¬ 
tal  F  factor,  but  tends  to  overestimate.  A  possible  explanation 
for  an  overestimated  F  factor  from  the  shear  map  is  die  differ¬ 
ence  between  the  altitude  of  the  aircraft  and  the  radar  beam. 
For  most  of  the  events,  the  aircraft  penetrated  die  microburst  at 
a  much  liiglier  altitude  dian  the  radar  beam.  Physical  observa¬ 
tions  .md  modclmg  results  suggest  diat  the  horizontal  shear  ui 
a  microburst  varies  with  altitude.  Tlius,  it  would  seem  pmdent 
to  attempt  to  compensate  for  the  discrcp.mcy  between  die 
hciglit  at  which  die  TDWR  antenna  beam  and  die  aircraft  mea¬ 
sured  die  microburst  intensity. 
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Figure  3  Shear  map  V3.  aircraft  horizontal  F  factor 
at  F  total  peak  time. 

Tlic  VicroyC  1991)  analyiical  microbursi  mode!  in¬ 
cludes  a  venical  shaping  function  for  die  horizontal  wind  ve¬ 
locity  diat  is  a  goodfittoexpcrimcntal  data.  Correctingforalti- 
tude.  Figure  4  shows  that  Uierc  is  a  marked  improvement  in 
die  shear  map  esdmates  for  horizontal  F  factor.  Therefore,  us¬ 
ing  the  shear  map  and  correcting  for  altitude  seems  to  provide 
an  acceptable  esdmation  of  the  horizontal  F  factor. 


6.  COMPARISON  OF  AIRCRAFT  AND  RADAR 

ESTIMATES  OF  VERTICAL  F  FACTOR 


From  formula  (1),  the  vertical  term  of  die  F  factor 
is  estimated  using  the  following  formula  (Again:  K’  is  equal 
to  one  because  the  shear  map  is  a  one  kilometer  shear): 


F 


=  AV/2[l\ 
ARlGSj 


(3) 


Figure  5  shows  die  shear  map  estimated  vertical 
term  versus  die  in  situ  F  factor.  A  probable  explanation  for  die 
poor  performance  of  die  downdraft  estimates  is  die  simplistic 
assumption  used  to  estj  mate  die  downdraft  velocity  Observa- 
uons  have  shown  dial  the  downdraft  velocity  varies  widi  alti¬ 
tude  as  well  as  across  die  radius  of  a  microburst.  A  better  es¬ 
timation  of  die  vertical  F  factor  needs  to  be  developed  diat  is 
capable  of  incorporatmg  the  aircraft's  locaUon  widiin  the  mi- 
croburst. 


Figure  4  TDWR  shear  map  vs  an  era  ft  horizontal  F  fac¬ 
tor  at  F  total  peak  lime  using  altitude  profile  correction. 


Figure  5.  TDWR  shear  map  vs  aircraft  i  eriical  F 
factor  at  F  total  peak  lime. 


7.  CONCLUSION 

While  die  current  TDWR  niicroburst  algorithm  per¬ 
forms  extremely  well  m  detecting  microburst  hazards,  some 
enhancements  are  needed  to  improve  its  ability  to  characterize 
the  hazard  m  terms  of  the  F  factor.  It  has  been  shown  that  die 
current  microburst  shapes  overestimate  the  F  factor  hazard  if 
die  aircraft  does  not  encounter  the  core  of  die  microburst.  A 
shear-based  approach  was  developed  which  allows  die  hori¬ 
zontal  F  factor  to  be  estimated  accurately.  However,  the  verti¬ 
cal  F  factor  term  remains  poorly  estimated  due  to  an  overly 
simplisdc  mass  continuity  assumption. 

In  orderto  improve  dieestimaie  of  the  vertical  F  fac¬ 
tor,  future  research  will  focus  on  fitting  an  analytical  micro- 
burst  model  to  shear-based  microburst  detections.  Such  a 
shear-based  algoridim  is  currently  under  development  as  part 
of  the  Integrated  Terminal  Weather  System  (TTWS)  program 
(Dascy,  1992).  Tliis  new  algoridim  will  allow  die  microbursi 
hazard  to  be  more  accurately  characterizcdbyprovidmgbctter 
localization  of  regions  of  intense  horizontal  shear  and  a  better 
esdmation  of  downdraft  intensity. 
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COHERENT  PROCESSING  ACROSS  MULTl-PRI  WAVEFORMS 


Mark  E.  Weber  and  Edward  S.  Ctaomoboy 


MIT/lJncoln  Lab 
Lexington,  Massachusetts 


1.  INTRODUCnON 

Meteorological  DoppIeriadarshave^callyutiUzed  constant 
pulse-rqietition  intervals  (PRI)  to  facilitate  clutter  filtering 
and  estinutdon  of  weather  echo  spectral  moments  via  pulse- 
pairorperiodogtam-basedalgorithms.  Utilixationofvariable 
PRIs  to  support  resolution  of  velocity  ambiguities  has  been 
discussed,  for  example  by  Banjanin  and  Zmic  [1],  but  not  im¬ 
plemented  owing  to  difficulties  associated  with  clutter  filter 
design.  Recent  work  by  Chomoboy  [2]  presents  de'^gn  algo¬ 
rithms  for  time-varying  finite  impulse  response  (FIR)  filters 
that  achieve  Chebyshev  or  mean-squared  error  (MSE)  op¬ 
timality  when  processing  multi-PRI  waveforms.  This  paper 
is  a  follow-on  to  duU  work,  treating  techniques  for  post-clut¬ 
ter  filter  processing  (e.g.periodogram  estimation)  that  areap- 
propriate  for  such  waveforms. 

Our  approach  involves  a  least-squares  fitting  of  the  signal — 
sampled  at  a  nonuniform  rate  —  to  a  weighted  sum  of  uni¬ 
formly  spaced  sinusoids.  The  sinusoids  or  “basis  fimctions” 
are  chosen  to  qian  a  Nyquist  interval  consistent  with  the  long¬ 
est  PRI  in  die  transmitted  waveform,  and  need  not  be  centered 
at  zero  Doppler.  Determination  of  the  sinusoid  weightings — 
effectively  adisaete  Fourier  transformation  (DFT) — and  the 
associated  residual  between  the  harmomc  fit  and  the  dau  are 
accomplished  viamuldplicationsofthesignal  vector  withpre- 
cotnpiued  ntatrices.  The  resulting  spectrum  estimate  can  be 
used  directly  for  weather  echo  moment  calculations,  or  can  be 
inverse-Foutier  transformed  usingconventional  techniques  to 
generate  a  time-domain  signal  representation. 

This  work  has  been  motivated  by  a  specific  application — es¬ 
timation  of  weather  spectrum  moments  for  a  Wind  Shear  Pro¬ 
cessor  (WSP)  modification  to  the  Federal  Aviation  Adminis¬ 
tration’s  Airport  Surveillance  Radar  (ASR-9)  [3].  Our 
approach  supports  candidate  low-altitixle  radial  wind  estima¬ 
tion  algorithms  [3]-{6]  thatoperateonfrequency-domainsig- 
lud  representations  and  require  that  the  radar's  block-stagger 
PRI  and  the  possibility  of  velocity  ambiguities  be  accounted 
for  in  generating  the  spectrum  estimates.  In  principle,  howev¬ 
er,  these  processing  techniques  are  also  applicable  to  weather 
radar  tystems  such  as  WSR-S8D  and  Terminal  Doppler 
Weather  Radar  (TDWR)  where  range  and  Doppler  ambiguities 
are  an  operational  concern. 

2.  LEAST  SQUARES  HARMONIC  FITTING 
Data  samples  at  arbitrary  times  t(n), 

x„=x[t(n)]  n  =  0.N-l  (1) 

are  modeled  as  the  weighted  sum  of  M  harmonically  related  si¬ 
nusoids; 
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ia=M-I 

=  2  ya,  exp{i2n(t)-HnAf?'.(n)j 

maO 

Weights  yg,  are  chosen  so  as  to  minimize  the  residual  between 
the  harmonic  fit  and  the  data  samples: 

£2  =  2:  I  Xa_x;|2  •  (3) 

iH) 

The  solution  and  conesponding  residual  are: 

y  =  ‘Sx  (4) 

£2  =  5H5x 

Hefey  is  theM-length  column  vectorofharmonicweightsyB, 
and  X  is  the  N-length  column  vector  of  data  samples.  The 
MxN  matrix  X  and  the  NxN  matrix  E  are  given  by 

**  (5) 

5=I-5IX 

where  basis  function  elements 

<j)n4n  =  exp  I  i2jt(fo+mAf)t(n)J  (6) 

define  the  NxM  matrix  C>. 

The  following  considerations  establish  the  number,  M,  and 
spacing,  Af,  of  sinsuoids  used  to  model  the  signal.  First,  the 
number  of  sinusoids  must  be  less  than  or  equal  to  the  number 
of  data  samples  so  that  the  system  of  equations  (2)  is  not  under¬ 
determined: 

M<N  .  (7) 

Second,  the  transform  of  the  frequency  sampling  “comb"  (a 
comb  in  the  time  domain  with  spacing  l/AQ  must  mt  fold  the 
signal  over  on  itself: 

Af<(t(NH(0)]-»  .  (8) 

Finally,  the  maximum  signal  bandwidth  representable  by  the 
sinusoid  set  must  be  consistent  with  the  longest  PRI  in  the 
transmitted  waveform: 

MAf  <  [max{t(n)-t(n-l))]"*  .  (9) 

We  choose  Af  so  as  to  satisfy  the  equality  in  equation  (8),  then 
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set  M  to  be  the  largest  integer  that  satisfies  equation  (9).  This 
guarantees  that  condition  (7)  is  satisfied. 

3.  APPUCAnONTOASR-9WSP 

The  ASR-9  utilizes  a  variable  PRI  to  mitigate  ‘'blind”  speeds 
for  aircraft  targets.  During  the  period  in  wluch  the  antenna 
scans  one  beamwidth  in  azimuth,  a  block  of  eight  pulses  is 
transmittedatalongPRI,foIIowedl7tenpulsesatashortPRL 
Because  the  antenna  rotation  rate  of  the  ASR-9  varies  under 
wind  loading,  “fDl  pulses”  at  the  long  PRI  may  be  inserted  fol¬ 
lowing  die  two  pulse  blocks  in  to  maintain  scan-to-ccan  azi¬ 
muth  re^stration  of  the  waveform.  The  ratio  of  the  long  and 
short  PRIs  is  9:7  with  a  Q^ical  value  for  the  long  PRI  of  1  ms. 
The  associated  Nyquist  interval  for  the  S-band  radar  is  53  m/s. 
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Figure  1.  ASR-9  transmitted  waveform. 

In  order  to  obtain  a  sufficient  number  of  samples  for  clutter 
suppression  and  Doppler  velocity  estimation,  tf.e  WSP  oper¬ 
ates  coherently  across  three  of  the  individual  pulse  blocks  as 
shown.  This  “extended”  coherent  processing  interval  (CPI) 
spans  27 successivepulsesandisthelongcstdeterminisiic  wa¬ 
veform  avulable,  owing  to  lack  of  a  priori  knowledge  as  to 
how  many  rill  pulses  will  be  insetted. 

Hgure  2  illustrates  a  csndidate  signal  processing  sequence 
used  by  the  WSP  to  generate  weather  moment  estimates.  In- 
phase  and  quadrature  signals  are  high-pass  filtered  using  the 
shift-variantFIRdesignsdescribedin[7].  Groupdelay  ischo- 
sen  so  that  the  filter  output  sample  spacing  is  equal  to  that  of 
the  input.  Two  samples  at  each  end  of  the  output  data  vector 
are  discarded  to  minimize  filter  degradation  at  the  beginning 
and  end  of  the  sequence.  A  rirstestimateofunambigousmean 
Doppler  velocity  is  obtained  through  application  of  the  “Chi¬ 
nese  Remainder  Theorem”  to  pulse-pair  Doppler  estimates 
obtained  individually  from  the  long-  and  short-PRI  data 
blocks. 


ter.  Weathermomentcsiimatcs,suchas“lowaltitude”Doppler 
velocity  [3]-{6].  are  generated  from  this  full-resolution  spec¬ 
trum  estimate. 


Figure  2.  Block  diagram  of  candidate  signal  processing  se¬ 
quence  exploiting  non-uniform  PRI  waveform  for  velocity 
amtnguity  removal. 


With  this  procedure,  unaliased  power  spectrum  estimates  are 
obtained  for  weather  signals  between  0.8  and  -0.8  times  the 
long-PRI  Nyquist  interval.  The  resulting  extended  Nyquist 
interval  (+/-  42  m/s)  is  sufficient  for  any  weather  conditions 
where  aircraft  landings  or  takeoffs  would  be  attempted. 

In  the  ASR-9  WSP  application,  resolution  cells  are  revisited 
onceeveiySseconds.  Since  weatherparameters  do  notevolve 
this  rajndly,  the  velocity  ambiguity  processing  need  not  be  re¬ 
peated  on  every  scan  of  the  antenna  (once  per  minute  is  sufri- 
cient).  As  shown  Ity  the  dashed  lines  in  Figure  2,  on  the  re¬ 
maining  scans  the  matrix  "K  can  be  preselected  and  cascaded 
with  the  matrix  TI  that  implements  the  shift-variant  FIR  rilter. 

Utilization  of  a  matrix  multiply  to  accomplish  the  signal  pro¬ 
cessing  leads  to  considerable  flexibility.  Forexample,  clutter 
riltering,  time-series  daia“windo  wing”  aitd  time-to-frequen- 
cy  trartsformation  can  be  achieved  through  a  single  matrix 
multiply  operation: 


y  =  e«WeiSH^  .  (10) 

Here,  the  elements  of  @  are  those  of  a  conventional  inverse 
DFT,  and  W  is  asquare”windowing”  matrix  whose  iron-diag¬ 
onal  elements  are  zero,  and  whose  diagonal  elements  are  the 
desired  wiixlow  function.  If  zero-padding  to  obtain  finer  ve¬ 
locity  spacing  of  spiectium  estimates  is  desired,  zeroes  are  ap¬ 
pended  to  the  columns  of  ©  and  the  rows  and  columns  of  W, 
and  is  replaced  tty  a  square  DFT  matrix  of  appropriate  or¬ 
der. 

4.  ILLUSTRATION  OF  MULTl-PRI  PROCESSING 


With  a  post-clutter  filter  data  vector  of  length  N=23,  the 
conditions  of  equations  (7)  through  (9)  establish  M=21  and 
Af = 48  s”*  (25  m/s).  Based  on  the  Chinese  Remainder  Theo¬ 
rem  estimate  c  ‘'unambiguous  mean  Doppler,  one  of  two  ma¬ 
trices  A ,  conesjronding  tobasisfunctionswith  center  frequen¬ 
cy  offsets  fo  of  respectively  plus  and  minus  03  times  the 
long-PRI  Nyquist  interval,  is  selected.  Additional  tests,  de¬ 
scribed  below,  may  be  used  to  confirm  that  the  selected  basis 
function  centerfrequency  offset  is  appropriate.  A“windo  wed” 
version  (see  below)  of  the  matrix  A  is  used  to  generate  a  spec¬ 
trum  estimate  from  the  full  23-sample  ouput  of  the  clutter  fil- 


Figure  3  shows  spiectrum  es  tim  ates  forasinsuoidof  frequency 
0.65  times  the  long-PRI  Nyquist  interval;  the  sinusoid  has 
been  sampled  at  27  points  corresponding  to  the  ASR-9  wave¬ 
form  illustrated  in  Figure  1 .  The  left  and  right  panels  are  esti¬ 
mates  gencr-sted  usingequation(10)withthe  above  basis  func¬ 
tion  offsets  fo  of  respectively  minus  and  plus  03.  The  filter 
matrix  H  has  been  chosen  to  be  all-pass  in  this  example.  The 
correct  choice  of  basis  function  results  in  a  spectrum  with  its 
peak  correctly  positioned  and  Doppler  sidclobes  that  are  con¬ 
sistent  with  theoretical  performance  of  the  Blackman  taper 
used  to  construct  the  window  matrix.  Incorrect  cl.oice  of  fo  re- 
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TROn,  OF  CORRECT  UNFOLD"  POWER  (dB) 


suits  in  signflcant  whitening  of  the  signal  spectnnn  and  an  or¬ 
der  of  magnitude  increase  in  the  residual  . 


DOPPLER  VELOCITY  (M/S) 


Figure  3.  Spectrum  estimates  &om  equation  (10)  applied  to 
simulated  sinusoid  sampled  atnoiv-equal  intervals.  Sinusoid 
normalized  firequency  is  0.65  and  basis  function  offsets  are 
-03  (left)  and  03  (right). 

Criteria  testing  the  appropriateness  of  the  choiceof  basis  func- 
tionftequencyoS'set  f^canbeusedtoreducethelikelihoodof 
a  gross  unfolding  error  from  the  initial  Chinese  Remainder 
theorem  Doppler  estimate.  As  illustrated  above,  examples  of 
such  criteriaarethemagnitudesoftheharmonicmodel  residu¬ 
als  and  the  “whiteness"  (power  ratio  of  minimum  to  maxi¬ 
mum  spectrum componeiu)  of  thespectra  estimated  using  dif¬ 
ferent  choices  for  f^.  Experiments  using  simulated  weather 
signals  with  varying  signal  to  noise,  signal  to  clutter  and  spec¬ 
trum  widths  have  indicated  that  the"  whiteness"  test  is  general¬ 
ly  more  robustc. 


Figure  4.  “Probability  of  correct  unfold”  as  described  in  the 
text  versus  signal  mean  Doppler.  The  leftpanel  treats  the  initial 
Chinese  remainder  theorem  estimate.  In  the  right  panel,  the 
"whiteness"  test  is  used  to  affirm  that  this  initial  choice  is  qj- 
propriate 

Figured  illustrates  the  concepL  Monte  Carlo  simulations  were 
run  using  a  weather  signal  of  moderate  spectrum  width  (5  m/s) 
and  low  signal-to-noise  ratio  (5  dB)  whose  normalized  mean 
Doppler  frequency  was  varied  from  0  to  -tO.8.  The  figure  plots 
the  probability  that  thepositive  basis  function  frequency  offset 
iscorrectly  selected.  (Note  that  forsignals  with  mean  Doppler 
less  than  0.2  minus  the  signal  bandwidth,  either  choice  for  f© 
is  appropropriate).  Our  initial  estimate  (solid  curve  in  the  left 


panel)  sets  thesignof  theselectedbasis  function  fi-equencyoff- 
setequallothatoftheChineseRemaindertheoremDoppleres- 
timale  and.  in  this  example,  is  incorrect  about  33%  of  the  time. 
This  estimate  is  chedted  against  the  choice  for  U,  that  mini¬ 
mizes  the  “whitest:*"  of  the  rcsulung  speemrnL  The  “white¬ 
ness"  test  compleme.’’  ^  the  initial  estimate  by  providing  near 
perfect  selection  of  basis  function  ft-equency  offset  as  long  as 
signal  Doppler  is  not  so  high  that  significant  power  is  outside 
theshiftedNyquistiriterval.  When  both  tests  are  ^lied  (right 
panel)  the  probabili^  of  selecting  the  apfnopriate  sign  of  f^ 
or  flagging  the  data  as  suspect  owing  to  disagreement,  is  near 
unityouttoanormalizeds)~nal&equeiKyof03.  At  higher 
signal  Doppler,  the  frequency  offest  selection  accuracy  de¬ 
grades  to  the  “baseline”  value  associated  with  the  Chinese  Re¬ 
mainder  theorem.  This  degradation  at  high  Doppler  magni¬ 
tude  could  be  eliminated  by  testing  of  additional  basis 
functions  with  larger  frequency  offsets  (e.g.  ±0.9). 

5.  OTHER  APPLICATIONS 

B  ase  data  degradation  produced  by  range  and  Doppler  ambi¬ 
guities  remain  a  fundamental  problem  for  weather  radar,  par¬ 
ticularly  with  systems  such  as  WSR-88D  and  TDWR  where 
automatedmeteorologicaldetectionalgorithmsareused.  The 
NEXRAD  Technical  Advisory  committee  recently  identified 
range/Doppler  ambiguities  as  the  second  highestpiiority  un¬ 
met  technical  need  (after  data  archiving).  Initial  Operational 
Test  and  Evaluation  (OT&E)  of  the  TDWR  in  Oklahoma  City 
has  likewiseillustrated  thatsecondtiip  weather  contamination 
and/oi  incorrectly  dealiased  radial  velocity  estimates  m$y  de¬ 
grade  the  operational  capabilities  of  the  radar. 

The  work  described  here  and  in  references  {2]  and  IT]  point  the 
way  to  processing  techniques  that  could  ameliorate  these 
problems  when  coupled  with  signal  waveform  changes.  Reli¬ 
able  resolution  of  Doppler  ambiguities  would  allow  foropera- 
tion  at  a  lower  average  pulse  repetition  frequency,  which  in 
turn,  wouldreducetheimpcctofrangefolding.  Whilethesig- 
nal  processing  requirements  arc  considerable  (approximately 
200  MFLOPS  in  our  ASR-9  application),  rapid  evolution  in 
digital  processing  hardwarecapability  makes  such  approaches 
feasible.  We  note,  for  example,  that  commercially  available 
single-board  arraypaocessing  cards  achievingffiisthrouglqiut 
are  available  for  under  520,000, 
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1  INTRODUCTION 

The  trade-off  of  range  vs.  velocity  ambiguity  is 
fundamental  and  operationally  significant  for  many 
S-  and  C-band  pulsed  Doppler  weather  radars. 
Transmission  schemes  using  multiple  pulse  repet'- 
tion  times  (PRTs)  [i.e.,  nonuniform  pulse  spacing] 
offer  the  potential  for  extending  the  unambigu¬ 
ous  measurement  range  by  resolving  intervals  of 
velocity  ambiguity.  Unfortunately,  multiple  PRT 
methods  can  be  problematic  with  low-elevation 
scanning  when  ground  clutter  removal  is  required. 
We  have  constructed  both  Chebyshev  and  mean- 
squared  error  (MSE)  design  algorithms  (Chornoboy, 
1993)  that  deal  with  design  in  the  complex  do¬ 
main;  the  MSE  algorithms  are  described  below. 

2  DESIGN  ISSUES 

Consider  the  design  of  a  finite  impulse  response 
filter.  For  an  N-coefficient  filter,  as  many  as  K 
designs  may  be  required,  where  K  is  the  number 
of  distinct  pulse  arrangements  of  length  N.  Since 
it  is  easier  to  consider  a  specific  example,  we  fo¬ 
cus  primarily  on  the  simple  case  shown  in  Fig.  1, 
that  of  an  alternating-PRT  scheme.  Here  two  fil¬ 
ters  are  “multiplexed”  in  the  sense  that  one  set  of 
coefficients^  Hi  =  [hio  •••  operates  on 

sequences  beginning  with  the  longer  pulse  interval 
Ti ,  and  a  second  set  (J/2)  operates  on  sequences 
beginning  with  Tj. 

Velocity  estimates  for  the  alternating-PRT  sig¬ 
nal  of  Fig.  1  can  he  obtained  by  using  the  single- 
lag  autocorrelation  method  known  as  Pulse  Pair. 
If  Ri  =  Ri^i)  represents  an  autocorrelation  esti¬ 
mate  obtained  from  pulses  separated  by  Ti,  then 
the  Doppler  velocity  (shift)  can  be  estimated  as 
Vi  =  -(A/4t7i)  arg(Ri)  ,  where  A  is  the  radar 
wavelength.  Similarly,  estimates  R2  and  V2  can 
be  obtained  corresponding  to  the  interval  r2. 

The  estimates  Vi  and  V2  “fold”  at  values  of  V 

*  The  work  described  has  been  sponsored  by  the  Federal 
Aviation  Administration.  The  U.S.  Government  assumes 
no  liability  for  its  contents  or  use  thereof. 

^Throughout,  superscript  “7^  is  used  to  represent  mu 
trix  trtmspose;  ,  complex  conjugate;  and  “t" ,  conjugate 
tr2msposc. 
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Figure  1:  Alternating-PRT  Filtering  Scheme. 
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Figure  2:  Ideal  Dealiasing  Using  AV  Method. 

determined  by  the  intervals  Tj  and  72.  If  Tj  ^  T2, 
velocity  folding  can  be  detected  and  corrected  by 
examining  the  phase  difference  between  Ri  and 
R2.  This  can  either  be  done  using  the  “difference” 
arg(Ri.^),  as  discussed  by  Zrnic  and  Mahapatra 
(1985),  or  the  difference  Vj  —  V21  as  considered  by 
Sirmans  et  al.  (1976).  The  latter  is  illustrated  in 
Fig.  2,  which  plots  ideal  relationships  for  the  case 
2ri  =  372.  The  upper  plot  shows  the  folding  pat¬ 
terns  for  Vi  and  V2  as  functions  of  true  velocity. 
The  lower  plot  illustrates  the  behavior  of  the  dif¬ 
ference  AV  =  Vi  -  V21  and  to  the  extent  that  this 
quantity  has  a  unique  mapping  to  intervals  of  true 
velocity,  ambiguities  can  be  resolved. 

It  is  desirable  to  have  a  magnitude  response 
over  the  extended  velocity  interval  that  is  free  of 
“blind  speeds”  so  that  measurements  Vi  and  V2 
v/ill  always  be  available.  However,  because  R  phase 
is  key  to  velocity  estimation  and  ambiguity  reso¬ 
lution,  it  too  is  an  important  aspect  of  the  fil¬ 
ter  design  process.  Figure  3  shows  a  staggered- 
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Figure  3:  Filter  Response  and  AV  Transfer  Func¬ 
tion  for  Staggered  Design  Without  Phase  Control. 


PRT  design  similar  to  that  achievable  by  adapting 
uniform  sampling  methods  [Banjanin  and  Zrnic 
(1991)].  For  an  assumed  3:2  stagger  spacing,  a  33- 
coefficient  filter  was  designed  to  provide  a  stop- 
band  half  width  equaling  0.04  V).  Although  a 
near  “flat”  magnitude  response  has  been  obtained 
(top),  there  are  intervals  where  the  phase  response 
deviates  significantly  from  linear  phase  (middle), 
and  the  effect  of  this  phase  error  on  the  ideal 
AV  transfer  function  (bottom)  is  near  devastat¬ 
ing.  Banjanin  and  Zrnic  (1991)  have  described  a 
r.iethod  that  would  work  around  those  areas  where 
the  AV  measure  is  most  impaired. 


and  S  is  defined  similarly  using  “sin”  instead  of 
“cos”.  Second,  let  at  represent  the  desired  output- 
sample  times  (i.e.,  group  delay)  for  the  filters,_and 
construct  output  vectors  Ct  and  St,  where  Ct  = 
[:osuo<rk  •  •  -coswAf-iffi]^,  etc.  Finally,group  “sin” 
and  “cos”  terms  together  to  form  complex  input 
'i  =  C-i-jS  and  output  ^t  =  Ct+jS^,  and  spec¬ 
ify  the  desired  magnitude  response  via  an  M  x  M 
diagonal  matrix  D.  The  ideal  output  response  for 
filter  k  is  D^t',  ^'’e  approximation  to  this  is  9Ht, 
and  the  approximation  error  is  f  =  SilHt  —  D^t- 
The  weighted  squared  error  S^QS  can  be  mini¬ 
mized  to  obtain  an  MSE  solution  for  Ht-  The 
M  'X.  M  real-valued  matrix  Q  is  used  to  introduce 
relative  weighting  for  pass-,  transition-,  and  stop- 
band  reg'ons.  We  have  found  it  useful  also  to  place 
a  constraint  on  the  maximum  gain  of  the  filter, 
which  can  be  included  via  a  term  such  as  HjHt- 
Although  the  above  provides  a  design  based  on 
minimizing  the  complex-domain  error,  it  may  not 
be  sufficient  because  it  does  not  permit  indepen¬ 
dent  control  of  phase  vs.  magnitude  error.  Let 
Ct  =  Et(w)  =  uat  —  ^"Hk  represent  the  phase  er¬ 
ror  for  filter  k.  For  |fi|  <  jr/2,  the  trigonometric 
inequality 


kt|  <  ^|sin(wo-i  -  ZWjt)!  (3) 

holds,  and  where  \'Ht  \  #  0, 

sin  ZTfib  =  Y^^kn  sinwr„  .  (4) 


3  MSE  DESIGN  EQUATIONS 


The  frequency  response  of  filter  Ht  is  defined  by 
N-l 

^  ,  (1) 

n=0 

where  T„  is  the  time  of  the  input  sample  (rel¬ 
ative  to  the  filter  start).  Let  Vt  =  Vt{u)  repre¬ 
sent  the  desired  output  response.  The  filter  de¬ 
sign  problem  is  to  find  coefficients  Ht  that  best  fit 
functions  e  to  the  desired  frequency  response 
Vt. 

It  is  straightforward  to  set  up  MSE  design  func¬ 
tionals  by  taking  M  discrete  frequency  samples  of 
Ht  and  Vt-  First,  define  M  x.  N  real-valued  ma¬ 
trices  C  and  5,  where  C  is  given  by 


C  = 


cos  Wo  To 


COSWoT;V_l 


(2) 


COSWaT-iTo  •••  COSWAf-lTAT-l 


Equations  3  and  4  can  be  combined  to  yield  the 
phase-error  constraint 


kl< 


mk\ 


sin uj{at  -  T„)  , 

n 


(5) 


which  can  be  used  to  force  |tt  j  small,  to  the  extent 
that  |7ft|  cooperates.  Define  the  Af  x  W  matrix 
et  = 


sinwo(cri  -  ro)  •••  siTwo((rt  -  ta^.i) 

sinww_i(cri  -  ro)  sinwAf_i(<ri  -  t/v_i)_ 

and  form  squared  error  term  ^fP^t,  where  = 
Ot-Ht  and  P  is  an  optional  weighting  matrix  for 
the  phase  error. 

An  error  functional  for  the  phase-controlled  de¬ 
sign  can  be  written 

f(Ht)  =  S^Q£  +  ^lP^t  +  H[Ht  .  (6) 
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Figure  4:  Filter  Response  and  AV'  Transfer  Func¬ 
tion  for  Staggered  Design  With  Excessive  Phase 
Control. 
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Figure  6:  Filter  Response  and  AV'  IVansfer  Func¬ 
tion  for  Block-Staggered  Design  With  Phase  Con¬ 
trol. 


Ti  Tj 

“•f  t  t  f  ttttt  I  t- 

I -  <ii - 1— >12— I 

Figure  5:  A  Block-Staggered  Sampling  Scheme 

This  is  quadratic  in  Hk  and  has  solution 
Ht  =  [3?e('i'tg<lr)  +  e'(pet  + 

=  [C'^QC  +  ^QS  eJ’PGi  + 

•  [C'^QDCk+S'^QDSk].  (7) 

Note  that  this  solution  only  requires  real-domaiu 
computation. 

At  the  extreme  P  =  0  (i.e.,  no  phase  control), 
the  design  of  Fig.  3  results.  If  P  is  instead  set 
very  large,  placing  high  priority  on  a  linear  phase 
response,  the  design  of  Fig.  4  results.  This  second 
design  approaches  the  “split  uniform  PRT”  filter 
discussed  by  Banjanin  and  Zrnic  (1991).  Linear 
phase  is  obtained  at  the  expense  of  introducing 
notches  (blind  speeds)  in  the  magnitude  response. 
For  the  alternating-PRT  signal,  magnitude  response 
blind/dim  speeds  vs.  nonlinear  phase  is  a  funda¬ 
mental  trade-off  which  no  optimal  design  can  fully 
overcome. 

An  additional  “design”  option  exists  however 
for  weather-radar  application  because  coherent  av¬ 
eraging  over  many  intervals  T)  (r2)  is  typically 
employed.  A  simple  extension  to  the  design  illus¬ 
trates  one  practical  potential.  Consider  the  block- 
staggered  signal  structure  illustrated  in  Fig.  5.  Pulse 
spacing  Ti  is  repeated  rji  times,  followed  by  in¬ 


terval  Tj  times;  after  which  the  pattern  re¬ 
peats.  This  signal  requires  K  =  fji  +  f}2  filter 
coefficient  sets.  As  with  the  above  alternating- 
PRT  example,  it  is  unlikely  that  exact  linear  phase 
can  be  achieved  without  some  compromise  in  mag¬ 
nitude  response.  However,  the  added  complex¬ 
ity  of  the  pulse  pattern  enables  an  improved  bal¬ 
ance  between  magnitude  and  phase  response.  Fur¬ 
thermore,  since  there  is  variety  in  the  filters  af¬ 
fecting  Ti  (22)  intervals  (across  the  confines  of 
one  block),  phase  and  magnitude  responses  can  be 
balanced  among  filters  by  “dithering”  (distribut¬ 
ing)  the  error.  Very  satisfactory  response  pro¬ 
files  can  result  as  shown  in  Fig.  6,  which  shows 
the  results  for  a  (4,4)  block-stagger  design  using  a 
phase-control  weighting  intermediate  to  that  used 
in  Figs.  3  and  4. 
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ANOMALOUS  PROPAGATION  ASSOCIATED  WITH  THUNDERSTORM  OUTFLOWS 


Mark  E.  Weber,  Melvin  L.  Stone  and  Joseph  A.  Cullen 
MIT/Uncoln  Lab  Lexington,  Massachusetts 


1.  INTRODUCTION 

Battan  [1]  noted  that  ducting  of  radar  energy  by  anom¬ 
alous  atmo^heric  refi:active  index  profiles  and  result¬ 
ing  abnormally  strong  ground  clutter  can  occur  during 
three  types  of  meteorological  circumstance:  0)  Litgc 
scale  boundary  layer  temperatue  inversions  and 
associated  sharp  decrease  in  moisture  with  height  — 
these  ate  often  created  by  nocturnal  radiative  cooling; 
(ii)  warm,  dry  air  moving  over  cooler  bodies  of  water, 
resulting  in  cooling  and  moistening  of  air  in  the  lowest 
levels;  Oil)  cool,  moist  outflows  hom  thunderclouds. 

In  contrast  to  the  first  two  types  of  anomalous  propaga¬ 
tion  (AP),  radar  ducting  associated  with  thunderstorm 
outflows  is  quite  dynamic  and  may  mimic  echoes  from 
precipitating  clouds  in  terms  of  spatial  scale  and  tem¬ 
poral  evolution.  While  non-coherent  weather  radars 
(c.g.  WSR-57)  are  obwously  suscqitible  to  false 
storm  indications  from  this  phenomenon,  Doppler  ra¬ 
dars  that  select  the  level  of  ground  clutter  supptesson 
based  on  "clear  day  maps"  may  also  fail  to  suppress  the 
AP-induced  ground  clutter  echoes.  Operational 
Doppler  radar  systems  known  to  be  suscqitible  to  tlus 
phenomena  are  the  National  Weather  Service’s 
WSR-88D  (Sirmans,  personal  communication)  and 
the  Federal  Aviation  Administration's  Airport  Surveil¬ 
lance  Radar  (ASR-9)  six-level  weather  channel  (2]. 

In  this  paper,  characteristics  of  thunderstorm  outflow¬ 
generated  AP  ate  documented  using  data  from  a 
testbed  ASR-9  operated  at  Orlando,  Florida.  The 
testbed  radar's  rapid  temporal  update  (4.8  seconds  per 
PPI  scan)  and  accurate  scan-to-scan  registration  of  ra¬ 
dar  resolution  cells  enabled  characterization  of  the 
spatial  and  temporal  evolution  of  the  AP-induced  clut¬ 
ter  echoes.  V/e  discuss  implications  of  these  pheno¬ 
menological  characteristics  on  operational  systems, 
specifically  the  ASR-9.  Algorithms  for  discrimina¬ 
tion  between  true  precipitation  echoes  and  AP-in- 
di  .ed  ground  clutter  are  discussed. 

2.  AP  MEASUREMENTS  WITH  ASR-9 

The  ASR-9  testbed  in  Orlando,  Florida  [3]  operated  at 
2.8  GHz,  transmitting  a  1  MW,  1  ps  uncoded  pulse  at 
an  average  PRF  of  1100  per  second.  The  antenna’s 
half-power  beamwidth  is  1.4*  in  azimuth  and  approxi- 

Thc  work  described  has  been  sponsored  by  the  Federal 
Aviation  Administration.  The  U.S.  Government  assumes 
no  liability  for  its  contents  or  use  thereof. 


mately  S*  in  elevation;  the  tilt  of  the  antenna  places 
the  peak  of  the  elevation  pattern  2*  above  the  horizon 
and  the  lower  half-power  point  on  the  horizon.  In- 
phase  and  quadrature  signals  are  sampled  at  intervals 
of  0.77S  ps  and  simultaneously  stored  on  high  density 
instrumentation  tape  and  processed  to  generate  real¬ 
time  displays  of  precipitation  refiectivity,  mean 
Doppler  and  spectrum  width.  Small-scale  divergent 
outflows  (microbursts)  and  outflow  boundaries  (gust 
fronts)  are  detec'^  by  automated  algorithms,  dis¬ 
played  and  tracked  in  real  time.  The  instrumented 
range  for  the  algorithms  extends  30  km  from  the  radar; 
maps  of  weather  spectrum  moments  are  generated  to  a 


range  of  111  km. 


Figure  1.  PH  scan  of  reflectivity  (quaiuized  in  NWS  "VIP" 
units)  during  an  AP-episo<L.  Range  rings  at  10  km  intervals. 


Data  from  six  separate  occurrences  of  thunderstorm- 
generated  AP  during  the  months  of  August  and  Sqi- 
tember  1991  and  1992  were  examined  for  this  paper. 
Significant  enhancement  of  ground  clutter  during 
these  episodes  was  observed  at  ranges  up  to  50  nmi  — 
the  strongest  returns  exceeded  65  dBz  equivalent  re¬ 
flectivity  factor  and  the  largest  AP  area  observed  was 
about  200  square  kilometers.  Duration  of  the  signifi¬ 
cant  AP-induced  ground  clutter  episodes  varied  from 
1.25  to  more  than  2.5  hours.  During  these  episodes,  in¬ 
dividual  "patches”  —  closed  regions  containing 
echoes  in  excess  of  35  dBz  equivalent  reflectivity  — 
varied  in  duration  from  a  few  minutes  to  the  lifespan  of 
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the  AP  qjisode.  Figure  1  shows  a  reflectivity  map  sistent  sources  ofstrong  echoes  with  repeatable  spatial 

generated  by  the  ASR-9  testbed  during  one  of  these  reflectivity  patterns. 

episodes.  Normal  ground  clutter  has  been  removed  by  Figure  3  compares  power  spectrum  estimates  from 

highpass  filters,  selected  using  a  "clear  day  map  [4];  AP-induced  ground  clutter  breakthrough  and  strati- 

the  remaining  echoes  are  precipitation  and  AP-in-  precipitation.  Spectra  of  the  AP-induced  echoes 

duced  ground  clutter ;  the  latter  echoes  ate  primarily  indistinguishable  from  normal  clutter,  consisting 

to  the  south  and  west  of  tiie  radar  at  ranges  greater  than  Qf  a  zero  mean  Gaussian  component  with  spectrum 

IS  km.  Note  that  some  of  the  AP  echoes  to  the  west  are  width  (0.7S  m/s)  consistent  with  antenna  scanmodu- 

contiguous  to  or  even  embedded  in  20  to  35  dBz  preci-  lation.  Precipitation  echoes  as  sensed  by  the  fan-beam 

pitation  echoes  to  the  west  of  the  radar.  ASR-9 — even  the  low  mean  Doppler  stratiform  rain 


Each  AP  episode  occurred  following  the  passage  of  a 
thunderstorm  outflow  boundaiy  (gust  front)  over  the 
radar  site.  In  the  above  example,  a  strong,  eastward- 
moving  gust  front  passed  over  the  ASR-9  about  30 
minutes  prior  to  the  dqiicted  scan.  Pre- and  post-gust 
front  rawindsonde  soundings  showed  significant  cool¬ 
ing  and  moistening  in  the  lowest  S(X)  meters;  maxi¬ 
mum  changes  in  temperature  and  dew  point  were  at  the 
surface  and  equalled  6  and  3  *C  respectively.  The  re¬ 
sulting  increase  in  radio  refractivity  gradient  below 
500  mis  shown  in  Figure  2.  TheAP  areas  are  observed 
only  in  the  sector  behind  the  gust  front,  implying  that 
the  superrefractive  environment  must  be  maintained 
along  the  entire  path  between  the  radar  and  the  ground 
scatterers  responsible  for  the  echoes. 

In  general,  the  patches  of  strong  AP-induced  ground 
clutter  appear  suddenly  (when  the  outflow  boundaiy 
has  passed  5-10  km  beyond  the  radar),  remain  approx¬ 
imately  constant  in  intensity  and  spatial  extent  for  a 
period  of  time,  then  dissipate  rapidly  over  the  entire 
affected  area.  A  characteristic  time  scale  fix  onset  or 
dissipation  of  individual  AP  patches  is  5  to  10  minutes. 
During  the  constant  phase  of  the  AP  episode,  echo  in¬ 
tensity  variation  in  time  is  small,  consistent  with 
scan-io-scan  fluctuations  in  normal  ground  clutter 
cross-  section  [4],  Specific  geographic  areas  are  con- 


DOPPLER  VELOerry  (M/S) 

Figure  3.  Power  spectrum  estimates  of  AP-induced  clutter 
(1^)  and  stratiform  precipitation  (right).  Dashed  lines  show 
theoretical  antenna  scan-modulation  spectrum. 

3.  AP-Induced  Ground  auttcr  Rejection 

Use  of  a  high-pass  ground  clutter  filter  in  all  range- 
azimuth  resolution  cells  would  eliminate  stationary 
clutter  breakthrough  caused  by  AP.  Such  filtering  may 
not  be  desirable,  however,  since  low-Doppler  power 
removed  by  the  filters  may  result  in  biases  in  weather 
reflectivity  or  mean  Doppler  estimates.  This  effea  is 
exacerbated  in  the  case  of  the  rapid-scanning  ASR-9, 
since  the  available  coherent  processing  intervals  are 
short  (8  or  10  pulses)  and  the  transition  bands 
associated  with  achievable  high-pass  filters  are  large 
14]. 

Tlie  ASR-9  weather  channel  and  the  WSR-88D  at¬ 
tempt  to  minimize  these  biases  by  utilizing  site-spe¬ 
cific  clear  day  maps  of  normal  ground  clutter  to  select 
the  minimum  level  of  clutter  suppression  necessary  to 
achieve  acceptable  weather  signal  to  clutter  ratios. 
With  tiiis  scheme,  a  "censoring”  function  should  be 
introduced  to  identify  stationary  ground  clutter  break¬ 
through  caused  by  abnormal  propagation  conditions. 
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Data  from  afliected  resolution  cells  are  flagged  and  can 
be  cither  disregarded  in  subsequent  processing,  or  re¬ 
processed  using  a  more  attenuating  high-pass  filter  so 
as  to  suppress  the  clutter  component  of  the  echo. 

We  have  examined  two  algorithms  for  discriminating 
between  AP  ground  clutter  breakthrough  and  actual 
precipitation  echoes;  both  depend  on  the  differing 
qjectral  characteristics  of  AP  and  precipitation 
echoes.  Direa  calculations  of  echo  spectrum  mo¬ 
ments  were  utilized  by  an  ASR-9  Wind  Shear  Plroccs- 
sor  [3]  to  cen-sor  AP-induced  clutter  breakthrough. 
The  censor  flag  was  set  for  range-azimuth  resoludon 
cells  where  the  mean  Doppler  velocity  was  less  than  1 
m/s  and  the  echo  spectrum  width  was  less  than  1.S  m/s. 
Spatial  consensus  filters  were  applied  to  the  censor 
flags  (e.g.  M-of-^  filters  along  the  range  axis  or  2-di¬ 
mensional  median  filters)  to  remove  "qteckle" 
associated  with  weather  moment  estimate  ernns,  par¬ 
ticularly  in  the  lower  intensity  AP  areas.  Figuredillus- 
trates  the  effect  of  this  censoring  process  using  the 
scan  shown  previously.  The  censoring  process  large¬ 
ly  removes  the  AP  without  significant  impact  on  the 
precipitation  echoes.  The  capability  to  remove  AP 
echoes  from  the  ASR-9 ’s  six-level  reflectivity  display 
was  favorably  received  by  the  Orlando  Air  TYaffic 
Control  team  during  operational  testing  of  the  Wind 
Shear  Processor  in  1991  and  1992. 


Figure  4.  As  in  Figure  1  but  with  AP-discrimant  enabled. 


An  alternate  approach  [2]  exploits  the  "inverse 
matched  filter”  characteristics  of  the  high-pass 
ground  clutter  filters.  AP-induced  clutter  echoes  will 
be  subject  to  large  attenuation  when  passed  through 
these  filters;  attenuation  of  weather  echoes  with  high¬ 
er  mean  Doppler  and  spectrum  width  is  much  smaller. 
A  power  threshold  lest  applied  to  the  ratio  of  clutter 


filter  input  and  ouqiut  can  effectively  discriminate  be¬ 
tween  AP-induced  clutter  breakthrough  and  precipita¬ 
tion.  Details  and  performance  examples  are  provided 
in  [2]. 

The  above  techniques  appear  sufficient  for  Air  Traffic 
Control  tq^lications  where  some  errors  in  the  exact 
intensity  and  areal  extent  of  precipitation  echoes  are 
tolerable.  Improved  performance,  useful  for  example 
in  hydrological  applications,  may  be  obtainable  by 
augmenting  these  single-gate  spectral  discriminants 
with  "expert  system”  knowledge  on  the  characteristics 
of  AP-induced  echoes  and  the  likelihood  of  supene- 
fraction.  Elements  of  such  a  system  would  include 
measurements  of  the  spatial  statistics  of  the  echoes, 
surface  temperature  and  humidity  measurements  — 
potentially  augmented  by  a  refractometer,  reliable  au¬ 
tomated  detection  of  outflow  boundaries  [S]  and 
knowledge  of  the  locations  of  ground  scatterers  likely 
to  be  illuminated  during  AP.  The  ate  specific  infonna- 
tion  necessary  for  this  last  element  may  be  obtained  tty 
means  of  detailed  terrain  maps  and  appropriate  ptopa- 
gatiort  models,  or  experimentally  through  accumula¬ 
tion  of  statistics  on  scattering  regions  from  many  AP 
episodes. 

4.  Summary 

Anomalous  propagation,  while  well  documented  since 
early  work  on  radar  meteorology,  remains  an  opera¬ 
tional  problem  —  particularly  when  it  occurs  in 
association  with  thunderstorm  activity.  This  paper 
discussed  spatial,  temporal  and  spectral  properties  of 
AP-induced  ground  echoes  associated  with  surface 
outflows  from  thunderclouds.  We  used  data  from  a 
testbed  ASR-9  to  demonstrate  a  reliable  spectral  dis- 
crimirumt  between  AP-induced  ground  echoes  and  ac¬ 
tual  precipitation  returns.  In  combination  with  addi¬ 
tional  sources  of  information  relevant  to  the  onset  of 
AP,  we  believe  that  performance  sufficient  to  support 
automated  hydrological  monitoring  in  the  presence  of 
diis  interferer  can  also  be  achieved. 
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1.  INTRODUCTION 

As  part  of  the  Aviation  Weather  Development  Program  of 
the  Federal  Aviation  Administration,  a  high  resolution  winds 
analysis  system  was  demonstrated  at  Orlando  International 
Aii}^  (MCO)  in  the  summer  of  1992.  The  purpose  of  this 
demonstration  was  to  illustrate  the  winds  analysis  capalrility 
possible  from  operational  sensors  in  the  mid  '90s.  An  impor¬ 
tant  part  of  the  design  of  this  system  was  the  development  of 
a  procedure  for  the  assimilation  of  Doppler  data  firom  multi¬ 
ple  radars.  This  procedure  had  to  be  able  to  automatically 
haixUe  regions  with  missing  data  from  one  or  more  radars,  as 
well  as  avoid  baseline  instability.  The  two  operational  radars 
scanning  the  analysis  region  were  the  National  Weather  Ser¬ 
vice  WSR-88D  (NEXRAD)  radar  located  approximately  65 
km  east  and  slightly  south  of  MCO,  and  the  MIT  prototype 
Terminal  Doppler  Weather  Radar  (TDWR)  located  7  km  due 
south  of  the  airport  The  base  data  from  these  two  Dopplerra- 
dats  wae  the  major  information  component  for  the  analysis 
system. 

Our  system  irtcludes  the  most  recent  improvements  in  the 
winds  analysis  portion  of  the  Local  Analysis  and  Prediction 
System  (LAPS)  developed  by  the  Forecast  Systems  Labora¬ 
tory  (McGinley  etal.,  1991).  LAPS  is  designed  to  run  locally 
on  systems  affordable  for  operational  weather  offices  and 
takes  advantage  of  all  sources  of  local  data  at  the  highest  pos¬ 
sible  resolution.  Our  implementation  for  the  airport  terminal 
region  is  called  the  Terminal-area  LAPS  (T-LAPS).  LAPS 
formerly  had  a  technique  for  the  assimilation  of  data  from  a 
single  Doppler  radar.  We  have  modified  that  technique  for  the 
assimilation  of  data  from  the  two  available  radars.  Our  ap¬ 
proach.  using  a  Multiple  Single  Doppler  Analysis  (MSDA) 
technique,  is  more  suited  for  unsupervised  operational  analy¬ 
sis  than  traditional  Dual  Doppler  Analysis  (DDA),  because  it 
is  able  to  handle  such  problems  as  incomplete  data  and  base¬ 
line  instability.  We  will  describe  the  T-LAPS  analysis,  with 
particular  attention  to  our  implementation  of  MSDA,  and 
give  some  examples  from  our  demonstration. 


*71ie  work  described  berewts  sponsored  by  tbe  Federd  Avution  Admiotsirt* 
tion.  The  Ihitted  Sutes  Government  cssusies  no  lubilrty  for  hi  cooteot  or  use 
thereoL 


Z  LAPS  WINDS  ANALYSIS  OVERVIEW 

The  LAPS  winds  analysis  (.Mbers,  1992)  uses  a  Bames 
(1964)  objective  analysis  scheme.  The  arudysis  acquires  a 
background  wind  field  and  recent  adnd  observations  in  the 
analysis  region,  and  produces  an  arutlyzed  wind  field  on  b3D 
grid.  LAPS  was  designed  to  be  computationally  efficient  and 
compatible  vnth  a  badeground  wind  field  provided  ty  a  pre¬ 
vious  arudysis  or  numerical  forecast  model.  For  the  demon- 
stradon  ordy  the  horizontal  winds  were  analyzed. 

The  stqis  in  the  analysis  process  are  as  follows; 

1 .  For  each  observation,  the  difference  between  the  u  compo¬ 
nent  of  die  observed  wind  and  the  u  component  of  the  back¬ 
ground  wind  at  the  grid  point  nearest  the  observation  is  com¬ 
puted.  Likewise,  a  difference  is  conqnited  for  the  v 
component. 

2.  At  each  analysis  point,  w.-ighted  means  of  the  u  and  v  dif¬ 
ference  values  are  computed,  t  form  a  correction  term  which 
is  an  estimate  of  the  vector  diffeience  between  the  actual  wind 
and  the  background  wind  at  that  point  The  wdghts  dqiend 
on  the  horizontal  and  vertical  distances  from  the  observation 
location  to  the  analysis  point  a  radius  of  influence  that  varies 
locally  depending  on  the  amUent  data  densi^,  and  sensor 
type. 

3.  The  correction  terms  are  added  to  the  background  wind  to 
form  the  analyzed  wind  field. 

Doppler  radars  measure  the  component  of  the  wind  only 
along  the  radar  beam.  Before  the  above  process  can  be  ap¬ 
plied  to  Doppler  radar  data,  the  Dopplerobservations  must  be 
transformed  into  vector  observations. 

3.  LAPS  SINGLE  DOPPLER  ANALYSIS 

This  section  details  die  process  by  which  Doppler  radar 
observations  are  brought  into  the  LAPS  arudysis.  The  idea  is 
to  transform  the  single  component  observations  fror..  a 
Doppler  radar  into  vector  quantities,  and  then  to  use  these 
vectors  as  additional  observations. 

The  steps  used  to  bring  Doppler  obsen’ations  into  the 
analysis  are  as  follows: 

1.  An  analyzed  wind  field  is  computed  using  the  background 
wind  field  and  the  non-radar  observations  as  discussed  in 
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Section  2. 

2.  The  wind  field  from  Step  1  is  adjusted  at  points  with  a 
Doppler  wind  speed  estimate.  At  these  points,  the  component 
of  the  wind  along  the  radar  beam  is  set  to  the  Doppler  value. 
The  perpendicular  component  is  undianged. 

3.  The  resulting  wind  vectors  at  points  with  a  Doppler  value 
are  considered  to  be  “radar  vector  observations”. 

4.  The  final  analysis  is  computed  from  the  origiiud  back¬ 
ground  wind  field,  the  true  vector  observations,  and  the  “ra¬ 
dar  vector  observations"  as  described  in  Section  2. 

4.  MULTIPLE  SINGLE  DOPPLER  ANALYSIS 

The  Multiple  Single  Doppler  Analysis  (MSDA)  devel¬ 
oped  forT-LAPS  is  a  simple  extension  of  the  standard  LAPS 
Doppler  analysis.  The  stqjs  used  to  bring  multiple  Doppler 
observations  into  the  analysis  are  as  follows; 

1 .  An  analyzed  wind  field  is  computed  using  the  background 
wind  field  and  the  non-radar  observations  as  discussed  in 
Section  2. 

2.  The  wind  field  from  Step  1  is  adjusted  at  points  with  a 
NEXRAD  Doppler  wind  speed  estimate.  At  these  points,  the 
component  of  the  wind  along  the  radar  beam  is  set  to  the 
Doppler  value.  The  perpendicular  component  is  unchanged. 

3.  The  wind  field  from  Step  2  is  adjusted  at  points  with  a 
TDWR  Doppler  wind  speed  estimate.  At  these  points,  the 
component  of  the  wind  along  the  radar  beam  is  set  to  the 
Doppler  value.  The  perpendicular  component  is  unchanged. 

4.  The  resulting  wind  vectors  at  points  with  at  least  one 
Doppler  value  arc  considered  to  be  “radar  vector  observa¬ 
tions". 

5.  The  final  analysis  is  computed  from  the  original  back¬ 
ground  wind  field,  the  true  veaor  observations,  and  the  “ra¬ 
dar  vector  observations”  as  described  in  Section  2. 

At  a  point  with  two  Doppler  wind  estimates,  the  mea¬ 
sured  radial  component  from  TDWR  will  equal  the  radial 
component  of  the  “radar  vector  observation”.  The  difference 
between  the  radial  component  measured  by  NEXRAD  and 
the  corresponding  radial  component  of  the  “radar  vector  ob¬ 
servation”  dependents  on  the  angle  between  the  two  radar 
beams.  When  the  angle  is  90°,  the  difference  is  zero.  As  the 
angle  decreases  to  0°,  the  difference  increases  to  the  differ¬ 
ence  between  the  the  TDWR  and  NEXRAD  measurements, 
and  at  0°,  a  “radar  vector  observation”  is  equal  to  the  single 
Doppler  “radar  vector  observation”  computed  from  only  the 
TDWR  data.  The  TDWR  data  were  chosen  to  follow  the 
NEXRAD  data  in  the  MSDA  process  since  the  TDWR  is  lo¬ 
cated  closer  to  the  Orlando  International  Airport. 

5.  DISCUSSION  OF  MSDA  AND  DDA 

In  traditional  dual  Doppler  analysis  (DDA),  a  wind  vec¬ 
tor  is  computed  at  each  analysis  point  with  two  Doppler  ob¬ 
servations.  The  resulting  wind  vector  exactly  agrees  with 
both  Doppler  values.  A'hen  the  two  radars  have  independent 
looks  at  the  wind  field,  defined  as  30°  or  more  between  the 
directions  of  the  beams,  DDA  generates  very  accurate  esti¬ 
mates  of  the  wind.  This  points  to  two  difficulties  that  arise 
with  DDA  in  an  analysis  system  which  must  produce  an  anal¬ 
ysis  at  each  grid  point.  First,  not  every  grid  point  will  have  a 


Doppler  return  from  each  radar.  Second,  when  the  two  radars 
do  not  have  indqxndem  looks  at  the  wind  field,  DDA  be¬ 
comes  numerically  unstable.  This  baseline  instability  gets 
progressively  worse  as  the  angle  between  the  radar  beams  de¬ 
creases.  The  first  difficulty  can  be  overcome,  but  will  result  in 
an  uxu-ease  in  complexity  relative  to  MSDA. 

MSDA  on  the  other  hand,  automatically  handles  incom¬ 
plete  Doppler  data,  and  does  not  have  a  baseline  instability. 
When  two  Doppler  values  are  available  at  points  whae  the 
two  radar  looks  are  ind^rendent,  the  “radar  vector  observa¬ 
tion”  is  very  close  to  the  wind  estimate  produced  by  DDA. 
When  the  two  radars  do  not  have  independent  looks,  MSDA 
produces  a  numerically  stable  “radar  vector  observation” 
with  one  high  quality  component,  a  Doppler  measurement 
The  other  component  is  derived  from  the  noit-Doppler  data 
sources.  At  points  with  only  one  Doppler  value,  MSDA  again 
produces  a  “radar  vector  observation”  with  one  high  quality 
component  Structural  constraints  imposed  during  the  im¬ 
plied  filtering  in  the  final  analysis  sVep  ensure  that  the  wind 
structure  in  each  of  the  sub-areas  blends  well. 

Our  implemenution  of  MSDA  was  developed  as  a  rapid 
prototype  for  this  demonstration.  As  such,  it  has  many  desir¬ 
able  properties.  However,  it  also  has  some  weaknesses  that 
we  will  address  in  the  future,  fri  regions  with  favorable  geom¬ 
etry  and  returns  from  both  radars,  the  “radttr  vector  observa¬ 
tions”  are  in  close  agreement  with  DDA,  but  are  then 
smoothed  by  the  aiutlysis.  When  the  two  radars  ate  looking  in 
nearly  the  same  direction,  the  NEXRAD  dau  are  largely 
overwritten  by  the  TDWR  data.  This  is  true,  for  example, 
even  when  the  analysis  point  is  closer  to  the  NEXRAD  than 
the  TDWR.  This  weakness  could  be  alleviated  with  a  weight¬ 
ing  between  the  two  radars  to  take  into  account  the  geometry 
of  the  analysis  region.  Each  “radar  vector  observation"  has  a 
different  level  of  quality  due  to  whether  the  observation  was 
built  from  one  or  two  Doppler  estimates,  and  the  radar  geom¬ 
etry  at  the  observation  location.  This  is  not  currently  taken 
into  account.  Lastly,  our  implementation  of  MSDA  can  be 
used  with  any  number  of  Doppler  radars,  but  even  with  3  or 
more  Doppler  radars  it  will  have  the  weaknesses  cited  above. 

Our  MSDA  implementation  is  equivalerd  to  producing 
“radar  vector  observations”  using  weighted  least  squares 
with  the  following  assumptions:  the  weights  are  inversely 
proportional  to  the  error  variance  of  the  data,  the  TDWR  error 
variance  is  infinitely  small  relative  to  the  NEXRAD  error 
variarree,  and  the  NEXRAD  error  variatKe  is  infinitely  small 
relative  to  the  error  variance  of  the  background  wind  field. 
This  suggests  optimization  as  a  path  to  improving  the  MSDA 
technique.  In  addition,  improvements  are  underway  to  gener¬ 
ate  a  background  wind  field  containing  pre-derived  dual 
Doppler  wind  vectors.  This  allows  the  existing  analysis 
scheme  to  improve  analyses  in  dual  Doppler  regions,  by  re¬ 
ducing  the  error  in  the  background  wind,  and  in  single 
Doppler  regions  by  increasing  tlie  accuracy  of  the  tangential 
components  of  the  "radar  vector  observations”. 

6.  TWO  EXAMPLES 

The  T-LAPS  analysis  region  was  120  km  x  120  km  in 
tlic  horizontal,  centered  on  the  Orlando  International  Airport, 
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and  extended  from  the  surface  to  a  height  of 500  nib.The  grid 
resolution  was  2  km  x  2  ion  x  SO  mb,  end  the  analysis  was 
performed  every  5  minutes. 

This  fine  grid  resolution  and  rapid  update  rate  were 
achieved  by  using  a  “cascadeof  scales”.  First,  the  winds  were 
analyzed  to  a  10  km  X 10  km  x  50  mb  grid,  every  30  minutes. 
All  of  the  available  data  sources  were  used  in  this  analysis, 
and  the  background  wind  field  was  derived  from  the  Mesos- 
cale  Analysis  and  Prediction  System  (MAPS)  (Benjamin  et 
al.,  1991).  Next,  the  final  analysis  was  performed  every  5 
minutes,  using  only  the  Doppler  data  and  automated  ground 
station  data,  LLWAS  and  ASOS/AWOS,  with  the  latest  10 
km  analysis  providing  the  background  wind  field. 

The  figures  show  the  analyzed  winds  at  400  fit  AGL,  and 
the  NfiXRAD  reflectivity  resampled  to  the  2  km  grid.  The 
wirxis  are  displayed  on  a  4  km  grid  to  reduce  visual  clutter, 
and  a  5  m/s  wind  arrow  is  shown  for  scale  in  die  upper  right 
comer  of  each  figure.  The  airport  runways  are  shown  in  the 
center.  The  four  outlines  are  lakes,  and  the  coast  appears 
along  the  northeast  in  each  figure.  Both  examples  are  from 
August  20, 1992. 

Figure  1  shows  the  wind  and  reflectivity  at  21:30  GMT. 
A  gust  front,  shown  both  a  reflectivity  thin  line  and  a  line 
of  convergence  in  the  wind  field,  is  being  produced  ty  a 
storm  ofi  the  coast  to  the  southeast  of  the  analysis  region.  Lat¬ 
er  in  the  day,  this  gust  front  collides  with  a  line  of  decaying 
storms  northwest  of  the  airport,  spawning  a  new  convective 
storm  system.  Figure  2  shows  the  wind  and  reflectivity  at 
23:55  GMT  associated  with  the  new  convective  storm. 


Figure  1.  Wind  and  Reflectivity 
(Aug.  20, 1992  21:30  GMT) 


Figure  2.  Wind  and  Reflectivity 
(Aug.  20, 1992  23:55  GMT) 


7.  ACKNOWLEDGEMENTS 

We  would  like  to  thank  the  Melbourne  National  Weather 
Service  Office  for  access  to  the  NEXRAD  base  data,  and 
Thomas  Schlatter  and  Stanley  Benjamin  at  NOAA/ERL/FSL 
for  providing  wind  data  from  the  Mesoscale  Analysis  and 
Prediction  System  for  the  background  wind  for  T-LAPS.  We 
would  like  to  thank  Russ  Bolton  and  Jerry  Starr  at  Computer 
Science-Raytheon  for  modifying  their  data  base  to  provide 
the  ground  station  data  in  our  analysis  region.  Finally,  we 
wculd  like  to  thank  the  LLWAS  program  office  for  expedit¬ 
ing  the  installation  of  the  LLWAS  3  anemometer  network  at 
the  Orlando  International  Airport 

REFERENCES: 

Albers,  S.C.,  1992:Thelaps  wind  analysis.  FourthWorkshop 
on  operational  meteorology,  Whistler,  B.C.,  Cana¬ 
da. 

Bames,  S.,  1964:  A  technique  for  maximizing  details  in  nu¬ 
merical  weather  map  analysis.  J.  Appl.  Met.,  3.  no 
4. 

Benjamin,  S.  G.,  er  a/.,  1991 :  An  isentropic  three-hourly  data 
assimilation  system  using  AC  ARS  aircraft  observa¬ 
tions.  Monthly  Weather  Review,  vol.  119,  no  4. 

McGinley,  J.,  S.  Albers,  P.  Stamus,  1991 :  Validationof  acom- 
positcconvective  index  as  defined  ly  areal -timelo- 
cal  analysis  system.  Wea.  Forecast.,  6,  337-356. 


39 


Adjoint-Method  Retrievals  of  Microburst  Winds  From  TDWR  Data* 
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1.  Introduction 

The  simple  adjoint  (SA)  method  of  Qiu  and  Xu  (1992, 
henceforth  referred  to  as  QX92)  was  recently  upgraded 
and  tested  with  the  Phoenix-II  data  for  retrieving  the  low- 
altitude  winds  from  single-Doppler  scans  (Xu  et  al, 
1993a,b,  henceforth  referred  to  as  XQY93a,b).  The  major 
results  can  be  briefly  reviewed  as  follows;  (i)  Using  mul¬ 
tiple  time-level  data  with  the  adjoint  formulation  makes 
the  retrieval  more  accurate  and  less  sensitive  to  the 
observational  error,  (ii)  Imposing  a  weak  nondivergence 
constraint  can  suppress  the  spurious  divergence  caused 
by  the  data  noise  and  improve  the  retrieval,  (ili)  Retriev¬ 
ing  the  eddy  coefficients  improves  the  wind  retrieval,  (iv) 
Retrieving  the  time-mean  residual  term  improves  the 
wind  retrieval. 

Although  the  results  in  XQY93a,b  were  encouraging, 
the  Phoenix-II  data  used  in  XQY93a,b  were  collected  on 
non-storm  days  with  chaff  dispensed  bom  an  aircraft.  The 
real  challenge  is  to  test  the  SA  method  with  storm  data. 
A  microburst  case  is  selected  for  the  test  in  this  paper. 

2. 11  July  1988  Microburst  case 

On  1 1  July,  a  very  strong  microburst  (>  35  m/s  differen¬ 
tial  velocity)  occurred  at  the  Denver  Airport  during  the 
1988  TDWR  (Terminal  Doppler  Weather  Radar)  opera¬ 
tional  test  and  evaluation  (Elmore  et  al.  1990,  Proctor 
and  Bowles  1992).  Dual  Doppler  coverage  was  provided 
by  the  TDWR  testbed  radar  (FL2,  operated  by  MIT  Lin¬ 
coln  Laboratory)  and  the  UND  (University  of  North 
Dakota)  radar  (see  Fig.  1).  The  operational  scan  strat¬ 
egy  executed  by  FL2  included  a  surface  sector  scan  over 
the  airport  every  minute.  This  surface  scan  was  matched 
nearly  simultaneously  (avg.  within  3.5  sec)  by  UND.  The 
polar  data  from  each  radar  were  threshold^  at  5  dB  SNR 
and  median  smoothed  with  a  5  gate  x  3  degree  filter  (at 
least  8  good  values  out  of  15  required).  The  data  were 
then  sampled  to  a  250  m  resolution  Cartesian  grid  (at  the 
level  of  z  =  190  m  above  the  FL2  radar  site). 

Surface  anemometer  data  from  the  12  station  Low 
Level  Wind  Shear  Alert  System  (LLWAS)  were  also 
collected  during  the  experiment  (see  Fig.  1).  Several  of 
the  stations  in  1988  suffered  from  wind  sheltering  prob¬ 
lems  (Liepins  et  al.  1990)  that  have  since  been  remedied 
by  raising  the  sensor  height. 


Fig.  I.  Locations  of  airport  runways,  radars  and  LLWAS 
stations.  The  inner  rectangular  domain  indicates  the 
region  where  the  winds  are  retrieved  in  Fig.  2a. 
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3.  Method  description 

As  in  XQY93b,  the  radial-component  wind  is  used 
as  a  "tracer"  field  and  is  governed  by  the  following 
approximate  radial-component  momentum  equation: 

9lVr  ♦  Vm-Vvr  -  v^m^/r  -  KVn^Vr  =  fm>  (D 

where  is  the  cross-beam  wind,  V  the  horizontal  vector 

wind,  (•)ni=(l/t^)Jo^(  )dt  the  time-mean  operator,  F 
the  unknown  residual  forcing  (mainly  the  pressure  gradi¬ 
ent  and  vertical  advection).  The  boundary  and  initial 
values  are  given  by  the  observed  Vf- 

The  objective  is  to  find  the  best  estimate  of  (Vm,  K, 
Fm)  in  (1)  that  gives  the  best  "prediction"  of  the  radial 
wind  Vf  in  terms  of  minimizing  the  following  cost-function 

J  E  1{P,A2  .  P2An,2  *  Pjdm^  -  P4U''»m-  (2) 

Here  {{(•))}E(l/Q)/J(  )dQ  is  the  area-mean  operator 
over  the  retrieval  domain  Q;  P|  and  P2  arc  nondimen- 
sional  weights,  AEV^-Vrob.  Am=Vrm"Vrobm-  and  ('^b 
the  observed  value  of  (•);  P3  and  P4  are  dimensional 
weights  (in  unit  m2),  djnsVH  Ym  *hc  divergence,  and 
the  vorticity.  The  minimum  of  J  can  be 
approached  by  numerical  iteration  along  the  gradient  of  J 
with  respect  to  (Vnj,  K,  Fm).  The  gradient  is  computed  at 
each  step  of  iteration  by  a  explicit  expression  derived 
from  the  adjoint  formulation  similar  to  (2.7)  of  XQY93b. 

The  optimal  retrieving  time  period  V  should  cover  4 
sequential  scans,  i.e.,  r=3Ar.  The  weights  are  given  by 
P,  =  lr/(t-At)r/2, 

Pj  =  0.02P|iii  with  Pim  =  (P|)m> 

P3  =  kjCJvr^Plm  "'‘d’  kj  =  30  ~  200  m2, 

P«}  =  k^Ovr^Pim  "’’dt  =  1(X)~  6(X)m2,  (3) 

where  Ovr  is  the  root  mean  square  amplitude  of  Vr.  The 
choice  of  the  time-dependent  form  for  Pj  was  explained  in 
QX92.  With  the  above  specified  value  for  P2,  the  weak 
form  of  the  constraint  A^i  =  0  can  reduce  the  error  in  the 
estimated  cross-beam  wind.  The  relative  strength  of  the 
weak  divergence  (or  vorticity)  constraint  is  controlled  by 
k3  (or  k4).  As  long  as  k3  (or  k^)  is  in  the  optimal  range 
shown  in  (3),  the  retrieval  is  not  very  sensitive  to  k3  (or 
k^).  The  weights  in  (3)  are  consistent  with  those  in 
XQY93a,b,  but  k/}  and  the  last  term  in  (2)  are  new  here. 

4.  Results 

The  SA  method  is  tested  with  the  microburst  data  for  a 
continuous  period  (22:04-22:33).  The  averaged  (over  25 
time-levels)  RMS  errors  and  correlation  coefficients 
between  the  retrieved  and  observed  variables  arc  listed 


*A  portion  of  this  work  was  sponsored  by  the  Federal 
Aviation  Administration.  The  views  expressed  are  those 
of  the  authors  and  do  not  reflect  the  official  policy  or  posi¬ 
tion  of  the  U.S.  Government. 
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in  Table  1.  When  the  observed  radial  winds  are  used  in 
the  final  results,  the  vector  RMS  errors  for  Vq,  reduce  to 
those  for  Vp<m  Table  1.  The  retrieved  wind  field  is  com¬ 
pared  with  the  observed  in  Fig.2a-b.  The  correlation  dia¬ 
gram  is  shown  in  Fig.  3,  where  the  RMS  error  and  corre¬ 
lation  coefficient  between  the  retrieved  and  observed 
wind  components  arc  also  listed.  The  retrievals  from  FL2 
radar  data  are  better  than  those  from  UND  radar  data. 

The  accuracy  of  the  retrievals  arc  affected  mainly  by 
three  factors:  the  data  noise,  the  temporal  fluctuation  of 
the  residual  forcing  (i.e.,  the  equation  error),  and  the  wind 
direction  relative  to  the  radar  beam. 

Using  the  wind  field  retrieved  at  the  previous  time  level 
as  an  initial  guess  can  reduce  the  CPU  cost,  but  may  not 
always  improve  the  accuracy.  Extrapolating  the  LLWAS 
data  to  the  grid  level  of  z  =  190  m  and  using  it  as  a  weak 
constraint  may  (or  may  not)  improve  the  retrieval,  if  the 
surface  winds  arc  well  (not  well)  correlated  to  the 
Doppler  radial  winds  at  the  the  grid  level. 


Table  1.  Statistics  of  the  retrievals  (with  FL2  radar) 


m/s 

IV 

'^(xm 

m/s 

^m  l^m 

10'3s- 1 IQ-^s- 1  1 0-2m/s2 

FL2  radar: 
RMS  error 
Correlation 

3.30 

0,92 

2.99 

0.83 

4.75 

0.60 

3.16 

0.22 

1.25 

0.77 

UND  radar: 
RMS  error 

4.53 

4.37 

5.34 

3.32 

1.41 

Correlation 

0,114 

..-0,0? 

0.48 

-0J.7— 

_ 

5.  Conclusion 

In  addition  to  the  earlier  findings  reviewed  in  section  1, 
it  is  found  in  this  paper  that  using  the  weak  vorticity 
constraint  also  improves  the  retrieval,  especially  for 
microburst  cases.  Using  the  previous  time-level  retrieval 
as  an  initial  guess  can  reduce  the  CPU  cost.  Optimal 
uses  of  the  surface  wind  data  need  further  investigations. 
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Fig.  2.  Comparison  between  the  (a)  retrieved  (from  FL2 
data)  and  (b)  dual-Doppler  observed  time-mean  wind 
fields  at  z  =  190  m  for  22:10-14, 16,  July  II,  1988. 
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Fig.  3.  Correlation  diagram  between  the  retrieved  and 
dual-Doppler  observed  winds  (for  every  5th  time-level 
during  the  period  of  22:04-33). 
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